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Abstract 

The problem under consideration has, actually, two aspects. The first one concerns the 
case when a nucleon is in free space that is in vacuum, while the second one studies a nucleon 
embedded into a nuclear environment. Both of these two aspects have been considered in 
the framework of chiral topological models of QCD. In this sense the whole content of the 
thesis may be divided into two main parts. 

In the first part the original Skyrme model (with m^r 7^ 0) has been extended by inclusion 
of the light scalar - isoscalar a - meson. The lagrangian has following advantages: 
It satisfies the scale invariance and trace anomaly constraint of QCD; 

It gives a realistic nucleon - nucleon (NN) interaction potential, which is quite close to the 
phenomenological one. 

The lagrangian has been further extended by explicit inclusion of a, p and uj - mesons as well. 
In order to get a more complete picture of NN potential the appropriate meson - nucleon 
vertex form - factors were obtained. 

The second part of the thesis considers a nucleon immersed into a nuclear medium. For 
this purpose a medium modified Skyrme lagrangian has been proposed. The lagrangian 
describes well such well known medium effects as decreasing of nucleon mass and increasing 
of its size ( swelling). We studied also a system with finite temperature (T 7^ 0) also. The 
temperature effects were taken into account by using the method of termofield dynamics 
(TFD). The explicit calculations show that the temperature dependence of meson - nucleon 
vertex form - factors is not simple. When the temperature increases they remain almost 
unchanged until a certain temperature Tc then start to decrease dramatically. The cor- 
responding critical temperatures for each meson - nucleon system were calculated and the 
appropriate conclusion about their origin was brought. 

The methods and results of this thesis may be used in the studies of heavy - ion collisions 
as well as of cosmology and of astrophysics. 
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Introduction 

The main goal of this thesis is the investigation of hadron properties in free space as 
well as in nuclear medium. It has been generally believed that, properties of hadrons and 
their interaction should be described in framework the of Quantum Chromodynamics (QCD) 
which is the underlining fundamental theory of hadron physics. However, in this theory one 
is failed with the difficulty related to the running coupling constant a(g^). In fact, for large 
momentum transfer g, as in the deep inelastic scattering of leptons by hadrons, a{(f) is small, 
and hence one can use the perturbation theory to make theoretical predictions which turned 
out to be in good agreement with the experiment. On the other hand, for the processes of 
interest in nuclear physics or low energy hadron physics, the length scale is typically ~ 1/m, 
that corresponds to small momentum transfer q. The running coupling constant a{q^) is then 
large and perturbation theory is useless. In this strong coupling regime, non perturbative 
methods are indispensable, but so far, not much successes has been achieved in this respect. 
In this regime one is led to modelling non - perturbative QCD. So, the challenge to nuclear 
physicist is thus to find models which can bridge the gap between the fundamental theory 
and our wealth of knowledge about low energy phenomenology. 

Our choice of models should be inspired by QCD, perhaps one day will be derived from it, 
and embody its important features, such as confinement and chiral symmetry, while avoiding 
the full complexity of a non - abelian guage theory. The parameters of these models must, 
for now at least, be determined phenomenologically, although the hope is that eventually 
they should be calculated from QCD. 

At present there is no systematic method for treating continuum QCD in the perturbative 
regime. With the exception of lattice Monte Carlo methods which still encounter both 



low - energy physics at the hadronic scale is mostly qualitative and relies heavily on the use 
of such phenomenological models as potential models [1], bag models [2], hybrid- bag models 
[3] and effective chiral models [4]. 

In the non relativistic potential models of Isgur and Karl [1] the quarks arc treated 
as constituent of protons much like in the conventional parton model. The confinement 
character of QCD is put in by hand using a linearly rising potential between the constituent 
quarks at large distances. The remaining gluonic effects are treated perturbatively, leading to 
Breit - Wigner hyperfine interactions. While these models are justified for the heavy fiavor 
hadrons such as charmonium and bottonium, they are questionable for the light hadrons 
such as N and A, where we know that, the constituent quarks arc highly relativistic. 

In the usual bag models, for instance, in the MIT bag model, the quarks are confined into 
a region of hadronic size V A^, inside of which the quarks carry their current masses and 
are treated fully relativistically. The bag is assumed to be a bubble of perturbative vacuum 
immersed into a non perturbative environment characterized by a bulk parameter B ( the 
bag pressure) [5]. Because of asymptotic freedom the remaining gluonic effects arc treated 
perturbatively inside the bag. Although rather successful in describing the spectroscopy of 
the hghtest hadrons, the MIT bag model suffers from a serious drawback, namely, axial 
- vector conservation. The confining bag wall breaks explicitly chiral symmetry, allowing 
left and right- handed quarks to mix at the boundary. We know that, in the chiral limit 
{rrif = ) QCD is chirally symmetric. The standard treatment of the MIT and cloudy bag 
model (CBM) lack translational invariance and they do not provide a dynamical description 
of hadron. This makes center of mass (cm.) and recoil corrections difficult to calculate 
[6, 7, 8, 9]. Besides, these models are not convenient to describe nucleon - nucleon dynamics 
as they make explicit use of the quark degrees of freedom. 

On the other hand, since the fundamental constituents in QCD, quarks and gluons are 
confined, it might be more suitable to eliminate their degrees of freedom which are not 
directly observable and incorporate their effects into an effective theory framed in terms 
of observable hadronic degrees of freedom including as many relevant features of QCD as 
possible. 

This alternative theoretical route can be traced back to the works of 't Hooft [10] and 
Witten [11]. 't Hooft proposed to look for the parameter in QCD, other than the coupling 
"constant" -a(g^) in terms of which one can perform expansions and arrived at the l/N^ 
expansion, Nc being the number of colors. Here one can note that, there is nothing funda- 
mental in the value Nc = 3; this is required by phenomenological considerations as the fit to 
the TTo 27 decay rate or to the ratio of e~^e~ cross sections, (T^Q^{e~^e~)/(T{e~^e~ fi'^^~). 

In a detailed analysis of QCD diagrams in the l/N^ expansion Witten [11] was led to 
the conjecture that in the large Nc hmit, QCD is equivalent to an effective field theory 
involving only mesons (and glueballs). Soon after this result, Witten [12] also recognized 
that baryons can be regarded as solitons of this meson theory in the manner first proposed 
by Skyrme [13]. 

All this results give support to the idea that low energy hadron physics can be described 
by effective theories constructed from meson fields alone. Therefore nonhnear chiral models 
including only meson degree of freedom will be used to study hadron properties in free space 
and in nuclear medium. Before entering into more details, we would like first to clearly 
identify what we mean by nuclear medium. There are obviously two different regimes, 
leading to two different physical representations of the medium. At very high density one 
expects a phase transition in which quarks and gluons are deconfined and more or less at the 
same time chiral symmetry is restored. Above the critical density Pc - which is expected to 
be of order of several normal nuclear matter density po - one thus expects a uniform medium 
of weakly interacting (almost) massless quarks and gluons. 



medium composed of more or less densely packed hadrons. Near the critical density it may 
be difficult to clearly identify nucleons in the medium dominated by all types of hadrons 
( mainly mesons and baryon resonances). At normal nuclear density or less the situation 
is however much simpler. Wc know that, the relevant degrees of freedom are nucleons 
interacting by the exchange of many correlated mesons [14]. 

We shall concentrate in this thesis on the low energy regime. One of the most exiting 
topics in this regime is the study of variation of hadron properties as the nuclear environment 
changes. Recent experiments from HELIOS-3 [15] and the CERES [16] at SPS/CERN 
energies have shown that there exists a large excess of e+e~ pairs in central S + Au collisions. 
Those experimental results may give a hint of some change of a hadron properties in nuclei. 
Ultrarelativistic heavy - ion experiments ( e.g. at RHIC) [17] are also expected to give 
significant information on the strong interaction (QCD) through the detection of changes in 
hadronic properties. 

Theoretically lattice QCD simulations may eventually give the most reliable information 
on the density and or temperature dependence of hadron properties in matter. However cur- 
rent simulations have been performed only for finite temperature (T 7^ 0). Therefore many 
authors have studied hadron masses in matter using effective theories and have reported 
that the mass decreases in nuclear medium. In the present thesis we shall make an attempt 
to modify Skyrme model to study the changes in mass of hadrons and in their interaction. 

The present thesis is organized as follows. In Chapter 1 the basic ideas and techniques of 
the Skyrme model "for beginners" is reviewed. In Chapters 2 - 4 we shall consider nucleon 
properties and nucleon - nucleon (NN) interaction in free space. Here the original Skyrme 
model is extended by inclusion of scalar - isoscalar meson which plays a crucial role in central 
the nucleon - nucleon interaction. 

It is believed that the nucleon properties and NN interaction may change in the nuclear 
medium. This modification is considered in Chapters 5-7. The main conclusion of the present 
thesis are summarized in the last chapter. 



Chapter 1 



The physics of Skyrme sohtons. 
(A brief review) 

An important property of QCD is its chiral symmetry. This invariance must be preserved 
in the corresponding effective field theories ^. 

1.1 The a model 

The simplest chirally invariant effective field theory is the a model proposed by Gell 
Mann and Levy [22]. Its version without fermions falls precisely in the class of effective 
field theories contemplated by Witten. This version of the a model can be formulated in the 
following way. Consider a 2 x 2 matrix field 

C/(r) = -^[(7(r)+if7f(r)] (1.1) 

where cr(r) is a scalar and isoscalar field, 7f(r) is a pion field and is the pion decay constant. 
The Lagrangian density can than be written as 

>C2(r) = \{d,af + \{d,i^f = ^Tr [d,U{r)d^U^{r)] (1-2) 

In the non-linear realization of the a model, the a field is eliminated from the model via 
the chirally invariant constraint, 

<^' + ^' = fl (1-3) 
implying that U is unitary. Defining a — fT^cos9{r) and 7? — fT^TtsinOlr), we can write 
f;_gir7r^(r) with tt = 7r/|7r| . 

The Lagrangian (1.2) possesses chiral SU{2)v ® SU{2)a symmetry. Left and right trans- 
formations are associated with left and right multiplication by elements of SU (2) group. The 
chiral SU{2)ji<Si SU{2)l group of transformations of the field U will correspond to the direct 
product of left and right transformations of the elements of SU (2) group: U — > AU with 
arbitrary constant SU{2) matrices A and B. The chirally invariant action is usually written 
by means of left (or right) invariant Cartan forms: 

^2 = -fTr(L,L^), 
L^ = U+d^U, R^ = Ud^U+. 
reader is referred to more extended review articles [18, 19, 20, 21] 



In terms of basic n and a fields the corresponding vector = [718^71] and axial currents 

— * 

Afj^ — adf^TT — irdf^a may be obtained by infinitesimal rotations in isospin space cr — > cr, 
TT — > TT — [en] and a ^ a — (en) , n ^ n + ea respectively. Clearly d^V^ — and d^A^^ — 
in accordance with chiral invariance. 

The trouble with the non-linear cr model is that it leads to energetically unstable solutions; 
by dimensional arguments, one can see that the energy scales like R, where R is the size of 
the system, and it vanishes when the system shrinks to zero size. 



1.2 The Skyrme Model 

To remedy this shortcoming, Skyrme [13] proposed the addition of higher-order terms in 
the derivative d^U to the form >C2(r). He suggested to introduce additional term of the form 

^4aW = ^Tr [{d,U)U\ {dM)U+r = ^Tr [L„ L^]\ (1.5) 

with the four derivative term. The resulting Lagrangian density 

^ C2{r) + C^^{r) (1.6) 

yields a complicate non-linear Euler-Lagrange equation which is difficult to solve. However, 
if one makes the ansatz that the pion field is directed radially in configuration space, n = r, 
f being a unit vector in coordinate space r = r/|r|, the U field is of the "hedgehog" form 
[3, 13] 

U = Uo = exp{irf9{r)). (1.7) 

In this case, the Lagrangian is 

L ^ j frC^, ^ ^ /; + ^^^^ 



with the choice of length scale x = 2efT^r and 6' = dO/dx. Minimizing L with respect to 9, 

dx y ae' ' ae 



i.e. using 4-{j^) = f| one gets the following Euler-Lagrange equation 



9 x9' sin 2^ 

^" — + 2sin2^ +- — + 

4 '24 (iq\ 

^,2 . sin2^sin2^ ^ ^ ^ 

+9'^sm2e = 0. 

x^ 

Soliton solutions can be found with the boundary conditions 9{oo) = and 9{0) = mr, n— 
integer. This solution is now energetically stable since the contribution of the term £4^^ of 
(1.6) to the energy is proportional to 1/R so that the total energy can be minimized with 
respect to R to give a finite value. It can also be shown [13] that the solutions characterized 
by this boundary conditions are topological solitons whose winding number is n. 



1.2. 1 The Baryonic Current. 

Besides the usual Noether currents, such as and A^, associated with various symme- 
tries fulfilled by the Lagrangian J0,g^{r), Skyrme showed that one can construct an "anoma- 
lous" current 

R ^ —}—p „_Tr r.'^rPr.'y (^ ^c\\ 



which is automatically conserved. A new quantum number can be therefore associated with 
Bf^ since this current does not correspond to any initial symmetry of the Lagrangian. With 
the hedgehog solution, the time component of is 

Integrating this density over all space r and taking into account the boundary conditions 

one gets B = J d^rBo (r) = n where n is winding number characterizing the solution U. At 

this point Skyrme identified n with the byron number, or equivalently B^ with the byronic 
current. Bciryons therefore emerge from the model as topological solitons. 



1.2. 2 Quantization. 

So far, the Lagrangian as well as the Euler-Lagrange equation are classical, t/ is a 
classical field and therefore carries no definite spin and isospin values. A simple quantization 
method starting from the classical solution U — Uq — exp (irri9(r)) was proposed by Adkins, 
Nappi and Wittcn [12]. These authors introduce collective coordinates through the unitary 
transformation A = 04 + ira where the independent of r but can depend on time t. 

The a^'s also satisfy the constraint a1 + (f = 1 since A is unitary. Under this transformation, 
Uo becomes U — A(t)UoA^{t). Substituting U into the Lagrangian yields 



L = y rfV>Csk[C/ = A{t)UoA+{t)] = -M + ATr {AA+) 



= -M + 2A^(a^) 



(1.12) 



where e.g. a = da/dt^ M is given by eq. (1.8) and 

A = — / dxx'' sin^ ^[1 + 4(e'2 + sin^ O/x^)] (1.13) 
3/7re'^ Jo 

Prom the Lagrangian one gets the Hamiltonian 

H = Il^d^-L, (1.14) 

where 's are the conjugate momenta = dhjddp, — 4Aa^ Substitution of this equation 
into the previous one yields 

^ = ^ + ^EnJ ■ (1-15) 

d 

With the usual canonical quantization procedure IIu = — ^7- — , the Hamiltonian takes the 
form 

4 Q2 

with the constraint a| + = 1. Because of this constraint the term 7—^ is to be inter- 

preted as the Laplacian on a 3-dimensional sphere. The wave functions of H are polynomials 
in a^. 



1.2. 3 The Spin and Isospin Operators 

Under a rotation about the z-axis in configuration space by a small angle S, it is easy to 
show that 

Hence, under this rotation the field configuration AUoA~^ transforms into a new one with 

(04 + ira) — > (04 + iTa)(l + iSTz/2) 

This gives after some algebra 

04 — > 04 — Saz/2, Qz dz + 5a^/2, 
a^:^ a^:- 5ay/2, ay ^ ay + 5aj:/2. 

Identifying this transformation with the initial rotation written in terms of the spin operators, 
i.e. — * (1 + i5Jz)a^ one finds 

% d d d d 

7, [-04^ h a^^ a^- h ay—] (1.20) 

2 oaz oa^ day oa^ 

Under a rotation about the z-axis of isospin space by a small angle S A ^ e' and 
the same calculation as before leads to 

id d d 

Iz^7,W4^ (1-21) 

2 oaz oa4 day oa^ 
It is easy to verify that the spin-isospin wave functions for various baryons are: 

b T) = (ox + «a2/)/7r , b i) = -^(a4 + «Oz)/7r , 

y 1 . 22 J 

\n t) = i{ai + ia^j-K , b i) = —{a^ — iay)/!: . 

The calculation of observables is then reduced to the evaluation of appropriate matrix el- 
ement of operators (constructed from the field configuration) with these spin-isospin wave 
functions. For example, the proton mass is 

(P }\H\p t) = ^ / d^a5{^al - l){a^-iay)H{a^ + iay) = 

o ^ ^ (1-23) 

where the expression (1.16) for H{U — AUqA'^) has been used. 



1.2. 4 Results for the Static Properties of Baryons 

The previous procedure has been appfied by Adkins, Nappi and Witten [12] to the calcu- 
lation of baryon masses (nucleon and delta resonance), charge radii, magnetic moments and 
coupling constants. The results are listed in Table 1.1. These authors adjust f.^ and e to fit 
the nucleon and A masses obtaining fj^ = 64.5 McV and e = 5.45. The other quantities are 
predicted. They are in qualitative agreement with experiment (within 30% ) except for qa- 



Table 1.1: Baryon properties in the original Skyrme model [12]. 



Quantity 


Prediction 


Experiment 


Mn 
Ma 


input 
input 
64.5 MeV 


939 MeV 
1232 MeV 
93 MeV 




0.59 fm 


0.72fm 


\' /M,I=0 

l^p 

fJ'n 

9A 

g-KNN 


0.92 fm 

1.87 
-1.31 
0.61 
8.9 


0.81 fm 

2.79 
-1.91 
1.23 
13.5 



1.3 Two baryon system 

The principal goal of studying the B = 2 sector of the Skyrme model is to apply it to 
the nucleon nucleon (NN) interaction. Characteristic features of the nuclear force shared by 
various semiphenomenological potentials are as follows: 

a) a long range (> 1.4 fm ) interaction dominated by one - pion exchange potential (OPEP), 

b) a medium - range (~ Ifm) attraction in the central force, and 

c) a strong short - range repulsion (< 0.7fm). 

The OPEP is well established, while the microscopic origins of the medium - and short range 
interactions are less well understood. 

A simple study has already been made by Skyrme himself, which surprisingly contains 
many essential features, at least qualitatively. He proved that a matrix product of two B = 1 
configurations, U = U1U2, belongs to the B = 2 sector. Then he proposed the "product 

— * 

ansatz" for two Skyrmions separated by R individually iso - rotated by A and B : 

U^^ = AUQ{f+ R/2)A+BUQ{f- R/2)B+ (1.24) 

The classical energies calculated using this ansatz giving 

E,2^2Eo + V{R,C), (1.25) 

where the static potential V depends of the relative distance R and the relative orientation 
C = A^B. A global rotation of U does not change the total energy. Skyrme pointed out the 
following two facts: 

i. For large R = \R\, V{R, C) approaches the one - pion - exchange potential (OPEP): 

V{R, C) ^ TabVaVb{eM-m.R)/R} (1-26) 

with Tab — Tr [TaCTbC'^]/2. The Yukawa function in eq. (1.26) is replaced by 1/R for 
the massless pion. This is natural since the Skyrmion contains a one pion cloud in the 
tail region. 

ii. At i? = and A — B — 1, the configuration U12 takes the form of hedgehog with B — 2 

and its energy is about 1 GeV higher than the two Skyrmion threshold. This excess of 



More quantitative analyses using the product ansatz have been performed by several 
authors and the precise form of the potential V{R, C) was worked out. The NN interaction 

— * 

is then obtained by the semiclassical method from V{R,C). Each Skyrmion is rotated by 
the collective coordinate A or B, and the relevant spin - isospin component is extracted. 
One finds that the NN potential consists of three terms, given by 

Vnn{R) = V^{R) + {m){B^a2)V,{R) + {m)S^2VT{R) (1.27) 

where 5'i2 is the standard tensor operator for nucleons: 

5i2 = 3(aii?)(a2i?) - (^1^2) (1.28) 

The last two terms in (1.27) tend asymptotically to the OPEP, as expected, while the first 
term has a shorter range. Unfortunately, it is found that the spin - isospin - independent 
part of the central force Vq is strongly repulsive at medium distances so that V/vat in Eq. 
(1.27) does not account for nuclear binding. 



d'S, = (^)F;,Fr, (1.29) 



1.4 Scaling behavior of the Skyrme model 

The second shortcoming of the Skyrme model was revealed by Shehter [23]. He proved 
that the chiral effective Lagrangian such as the Skyrme model have different scaling behavior 
than QCD. Classically the QCD Lagrangian (with massless quarks) is invariant under scale 
transformations with A real. The corresponding conserved current is the 

dilaton current, — Ti^^x^ where T^^ is the energy-momentum tensor. At the quantum 
level, however Sj^ has nonvanishing divergence [24] 

^/%). 
2^ 

where g is the coupling constant, I3{g) the QCD beta function and F^^^ the gluonic field 
strength tensor. Since d^S^ = the trace of T^^^ in Eq. (1.29) is called the "trace anomaly" . 
To restore the anomalous scaling behavior of QCD at the effective Lagrangian level one 
introduces a scalar field x(a;) [23] such that d'^S/^ — —bx^- Fluctuations in x around Xo niay 
be interpreted as a scalar glueball or a quarkonium. 

In the next chapter we shall take into account scalar dilaton field to restore the anomalous 
scaling behavior and consider NN interaction. Meanwhile, to make this point clearer and for 
further references we briefly outline the meaning of "scale dimension" below. Let's consider 
the kinetic terms of a fermion field (e.g. quark field in QCD) 

4in(^) = (1-30) 

and a boson (pion) field 

>^kinW = Id.nd'^n. (1.31) 

The scale invariance demands that the action S — J d'^xC should be invariant under scale 
transformations Xn — e^'^x^ : 

S^J d'xC]^^{x) = J dV4i^(xO (1.32) 

Let the fermion and boson fields scale as ip{x) — exp {nX)ip'{x') and Tr{x) — exp (mA)7r'(a;') 
respectively. For fermions we get 

,A„2nA„-4A _ 

(1.33) 



S^iJ d^xV'(x)9^7^V'(^) = ^ / dVV''(x')9'^7^V'(^')e^e"'""e 



and hence — 4A + A + 2nX = A(2n — 3) = 0. Since A is arbitrary then n = 3/2. Similarly, for 
boson fields one can show that m — 1. Therefore we may conclude that 



^x^, ipix) e^^/'^ip{x), 7r{x) e^7r{x). (1.34) 



i.e. the scale dimension d equals to df — 3/2 and db — 1 for fermions and bosons respectively. 
Obviously, the scale dimension of chiral nonlinear field U{x) is zero. For this reason the four 
derivative term of the Skyrme Lagrangian £4^^ is scale invariant by itself: 

j (i^x/:4a(x) ^ j d''x'e~^^C'^^{x')e'^^ = j d^x'C'^^{x') (1.35) 



whereas the kinetic term 



C^ix) = ^Tr [d^U{x)d^'U+{x)] (1.36) 



is not: 



£2 

I d^xC2{x) -^ -f j d^x'e-^^e^\Y[d^U\x')d^'U'+{x')] 
j d^x'C'^ix') 



(1.37) 



How can the invariance be restored? According to Shehter [23] an additional boson field 
x(x) should be included into the Skyrme Lagrangian. As a result the minimal version of 
scale invariant Skyrme Lagrangian will be given by 

2 4 xo (1.38) 

where xo is the vacuum expectation value of a dilaton field x{^)- 

But unfortunately this is not the whole story. The basic QCD Lagrangian whose main 
features should be reflected in an effective one has the "trace anomaly" mentioned above. 
To restore the anomalous scahng behavior of QCD at the effective Lagrangian level one has 
to introduce an effective potential V{x)^ describing the self interaction of dilaton fields. It 
was shown that, the effective potential could be uniquely determined by the vacuum energy 
density ^Xo/"^ condition that the minimum should occur at x = Xo which yield [24] 

V{x) = ^{1 + (^)l41n(x/xo) - 1]}. (1.39) 
The parameters h and xo are obtained from the value of the gluon condensate Cg 

mS,\^) = (0|(^)F;,Fr|0)^Cg = 6x^ (1.40) 



and the dilaton mass 

K=g^\x=xo-^bxl (1.41) 

which together with /^r and e specify the parameters of the Lagrangian. It has been shown 
that, the suitable choice of parameters Xo ~ l20MeV, Cg ?a (300Mey)^, e ^ 4.5 and 
/tt = 93MeV leads to a large dilaton mass ^ ISOOMeV", which may be identified with a 
glueball. Although this scalar field has nothing to do with a scalar a meson {m^ ^ 600Mey), 
its inclusion into the Skyrme Lagrangian proved to solve a long standing problem of missing 
of intermediate ra,np"e attraction in the central NN interaction f25. 2fil. 



Nevertheless, it would be quite desirable to include into the Skyrme Lagrangian a dilaton 
field with lower mass e.g. — rria ~ 600MeV. This goal has been achieved by using the 
method of bosonization in ref.s [21, 27]. As a result, a modified Lagrangian which differs 
from Shehters one [23] mainly due to the self interaction potential: 




(1.42) 



has been proposed. We shall use this model to study A^A^ interaction in the next chapter. 



1.5 PC AC and Goldberger - Treiman relations 

The chiral invariance of Skyrme model leads to well known relations. Here we show how 
these relations can be derived in a rather model independent way. Considering pion - nucleon 
system we start with the interaction Hamiltonian: 

Hint = -i9nNN^ix)'jr^f^{x)(^^{x) (1.43) 

where Qt^nn is the pion - nucleon coupling constant, ip{x) and (fini^) are nucleon and pion 
fields respectively. The latter satisfies Klein - Gordon equation: 

(d^d^ + ml)ifi^(x) = -ig^NNi^{x)l5ril^{x) (1.44) 
Axial current of this system is the sum of the nucleonic, ^4^, 

and pionic, A^^, 

a; = Ud^ifA^) (1.46) 

axial currents: 

A,^A^ + a; (1.47) 
Using Dirac {i^nd^ — M)'4>{x) — and Klein - Gordon equations one may obtain 



df^Al = -Uml0^{x) 



(1.48) 



Note that, the last second equation in (1.48) is called as PCAC (partial conservation of axial 

— * 

current) relation. Now, the requirement of conservation of the total axial current d'^A^ — 
gives following equation for the divergent of A^: 

{d^,d^ + ml)(p^{x) = -gAMi,{x)i-i^Tij{x) (1.49) 

Finally, comparing right hand sides of the Eq.s (1.44) and (1.49) one easily gets Goldberger 
- Treiman relation: 

MgA = UgnNN (1-50) 

This utmost important relation between axial coupling constant gA and pion - nucleon 
coupling constant gnNN is in good agreement with the experiment: g-„NN = 13.4 ± 0.1, 
MgA/ fir — 12.78 ± 0.1. The small deviation refiects the accuracy of the relation. We shall 
come back to this point considering in - medium modification of this relation in Chapter 5. 



Chapter 2 



Scalar Dilaton - Quarkonium Meson in Nucleon - 

Nucleon Interaction. ^ 

2.1 Introduction 

Skyrme - like models allow for the study of the properties of individual baryons and 

offer a framework for investigating interactions between classical baryons. The simplicity of 
calculations and the existence of a transparent connection between NN and NN interactions 
within the Skyrme model are distinct advantage over the models such as nonrelativistic quark 
models or the bag models. The most of calculations based on the product approximation 
within Skyrme - like models lead to projected potentials that are in qualitative agreement 
with empirically motivated NN potentials at short and long distances. However many of 
them does not take into account the scale anomaly of the QCD omitting scalar particles and 
fail to reproduce the intermediate range attraction in the scalar channel. 

At the present time there are two different approaches to include scalar meson - dilaton 
into the Skyrme model. In both approaches the interaction of dilaton field with the chiral 
field is dictated by the scale invariance. The main difference is in the origin of dilaton. On the 
one hand, there is the approach in which dilaton is associated with the glueball [31]. In this 
approach the glueball field saturates the scale anomaly completely. In the other approach 
the dilaton is treated as a quarkonium arising due to fiuctuations of the quark condensate 
[27], which is an order parameter for the chiral symmetry breaking. In this approach the 
dilaton - quarkonium saturates a part of the scale anomaly, namely by the part which is 
produced by the quark loops. This choice of the quarkonium as a dilaton is based on two 
principal observations: 

a) The experimental studies of scalar resonances [32] show that the real lightest candidate 

for the glueball state is /o(1590) which does not appear in tttt and KK productions. 

b) Consideration of the chiral anomaly [27] shows that the only guage - invariant combi- 

nation of the gluon field G^^Gxp can contact only with a total antisymmetric tensor 
^fiu\p^ Thus this combination has J^'" = 0^^ quantum number and contributes to 
Wess - Zumino - Witten action producing U{1) anomaly. It suggests that the O"*""*" 
glueball field as a fiuctuation of cannot interact directly with chiral field, but only 
through the mixing with the dilaton - quarkonium. The estimate gives small values 
for the mixing angle of glueball and quarkonium states. 

In the present chapter we calculate the central NN potential within a Skyrme model 
suggested by Andrianov and Novozhilov in ref. [27]. The model is based on effective low 

^The present chapter is based on following articles by the author and his collaborators: [28, 29, 30] 



energy action for pseudoscalar and scalar (dilaton - quarkonium) meson which has been de- 
rived starting directly from the QCD generating functional by the joint chiral and conformal 
bozonization method. 



(2.1) 



2.2 The Lagrangian 

The Lagrangian suggested by Andrianov and Novozhilov [27] is 

jC2^-^[TrL,L'^-2{d,aY]e-'^, 
Aa = ^Tr[L„L1^ 

W^. = -^[e-4^ -1 + ^(1- e-^^)] 

where /^r = QSMeV is pion decay constant, e— dimensionless parameter, Cg - gluon con- 
densate parameter : Cg = ((300 — 400)Mey)^, a - scalar meson field , — U'^d^U , U- 
is an SU{2) matrix chiral field. It is a generalization of the well - known original Skyrme 
model [12, 13] (OSM) and takes into account the conformal anomaly of the QCD. The term 
C2 includes the kinetic term of the chiral field and scalar fields. The C/^g^ is a well known 
Skyrme term. The effective potential for the scalar field is a result of the extrapolation of 
the low energy potential into the high energy region in one - loop approximation to the Cell 
- Mann - Low QCD (3 - function. The parameter e depends on the number of flavors Nf: 

e = 8Nf/[33 — ^Nf). Note that in the limit of heavy a- meson the potential becomes equal 

1 2 

to the symmetric quartic term = -^-^[Tr Lf^L^] , which is necessary to reproduce the tttt 

scattering data [33]. We extend the Lagrangian (2.1) by adding chiral symmetry breaking 
term: 

£^3b = 7^^T r {U - 1) (2.2) 
to provide the correct asymptotic behavior of the soliton field. Thus our model is given by: 

^ = ^AN + ^xsh ^2-2) 
In the static approximation when some parts of the Lagrangian have simple form 
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(2.4) 



where p = exp (— (T(r)) wc may use the most popular hedgehog ansatz: U = Uq = 
exp {iff6{r)), f = f/\r\ and propose spherical symmetry for the scalar field (j{r). Then 



the mass functional for classical soliton is given by: 

Mh = ^ dx{M2 + M/^^ + Mw + Mtt), 

M4a = s^(d/2 + n, (2-5) 

where x = 2e/^r, c = cos(^), s = sin(^), d = {s/x)'^, f3 = m7r/(2/^e), D^q — Cg/384e^/^^. 
The Euler - Lagrange equations for the chiral angle 9{x) and scalar meson shape function 
p{x) follow to be 

9"[^x^p^ + 2s^] + i [2x^9' pp' + 2x9' p^ - p^S2]+ 
+S2{9'^ -d) - p^x^Ps /4 = 0, 
x^p"/4: - px'^{9'^/2 + ci)/2 + xp'/2 - W^Qx'^{p^- 

-p^ - 1) - 3p2(l - c)x2/3V4 = , ^^'^^ 

where S2=sin(2^) and a prime corresponds to the derivative with respect to x. For the 
asymptotic behavior of 9{x) at large distances when p — > 1 we get from (2.6) and (2.7) a 
familiar formula: [34] 

e{x) ~ + exp (-/3x)^^^l±^ exp {-(5x) (2.8) 

where is the nucleon mass, Qttnn = 13.5 - the pion-nucleon coupling constant. This 
formula gives connection between gj^NN and Skyrme parameter e and reveals that e is not a 
free parameter of Skyrme model in some sense. The behavior of p{x) at large x exhibits a 
rapid decrease from the unity: 

—f)x 

p(^) ^ 1 _ , D ^ ^l6L>gff(4 - e) (2.9) 

At small distances the behavior of 9{x) and p{x) is given by: 

9{x)^7:-ax 
p ~ Po ~ ox 

with a = [3q;Vo + i^D^siPo^ - Po^ ~ ^) + Wf^V^- 

Thus Eq.s (2.6) and (2.7) may be solved numerically using well known methods [35. 36] 
and the parameters A, a, 7 and po in Eqs. (2.6), (2.10) may be adjusted to yield continuity 
in 9,p,9',p' at a merging point. In addition, the consequence of the virial theorem: 



f 

Jo 



dx[M2 - M4a + 3{Mw + Mn)] ^ (2.11) 



may be used to control the accuracy of stable solutions of Eq.s (2.6), (2.7). The resulting 
solution, 9{r), is similar to that one in the original Skyrme model (OSM) [12] while p(r) 
rapidly goes to the unity starting from Pq at r = . 

Masses of nucleon Mjv and A - isobar - Ma may be calculated from the following expres- 
sions: 

Mn^Mh + 3/8Am ^2 12) 



Ma^ Mh + 15/8 Am 
where Mh - soliton mass (5), Am is momentum inertia of the rotating skyrmion: 

47r 



Am = ^/ dxxV[p74 + ^'2 + d] (2.13) 



2.3 The nucleon — nucleon interaction 

We investigate the skyrmion-skyrmion interaction by using the product approximation: 

U{x]ri,Ai,r2,A2) = AiUo{x-fi)AtA2Uo{x-r2)Ai = U1U2 
p(f , rl, rl) = p(f - rl)p(f - rl) = P1P2 

where f/o(x — f^i) for i = 1,2 is the hedgehog solution (?7o(^) = exp (irr6'(r)), f = r/\r\) 
located at f^, and Ai the collective coordinate to describe the rotation. The skyrmion- 
skyrmion potential is defined by: 

V{r)^-j dx[jC{U,U2, P1P2) - Pi) - ^2, P2)] (2.15) 

where r is the relative coordinate between two skyrmions (r = — r^). The product 
approximation has reasonable properties: 

a) the baryon number is additive, 

b) it gives a correct behavior at the asymptotic region, i. e., the one-pion exchange potential, 

c) at r = the product form is very close to the hedgehog solution of the baryon number 

5 = 2. 

The static NN potential may be obtained by using a standard technique [37, 38, 39, 40] 
which gives the following expressions for the scalar-isoscalar part of central NN interaction: 

Vcir) = ~~^J^^ dxx'^[V2{x, xi) + V^aix, Xi) + Vw{x, Xi)+ 

+V^^^{x,xi) + 1 ^ 2] 

Here the contributions from C2 > ^Aa. ' '^xsb denoted by V2, V40, V^sb and Vw 

respectively: 

V2{x,Xi) = {p{x)(f){x)pi(f)iZi + (pi^ - l)x 

x[p\x)(e{x) + 3d{x)) + (f>\x)]}/8 

Vw{x,x,) = D^q{[p\x){p^^ - 2) + 1] - ^p%x){p,' - 2) + l]/£}/2 (2.17) 
V^a{x,xi) = {di[3d{x) + 29{x)] + 9{x)9i{l - ^i^)/2}/3 
V^^-^{x,xi) = -p^{x)d\x){pi^[c{x)ci - l]/2 - c{x) + l}/4 

where the following notations are introduced: ^^(x') = 9''^{x) — d{x) , (j){x) = p'{x) , 
(2^1,2)^ — x^ + r^±2xtr , Zi^2 — {x ± ft) /xi^2 , f — '^f^f^n , r is a relative distance between 
the two nucleons, pi = p{xi), di=d{xi), {i — 1,2) e.t.c. 



2.4 Results and discussions. 

As a starting point we make an attempt to describe the pion - nucleon coupling constant 
QttNN and N — A mass splitting: Mjva = — optimizing the only free parameter of 
the model - Cg. In fact, the value of /^r should be fixed to its experimental value, and the 
Skyrme parameter e must be defined by Eqs. (2.6) - (2.8) with Qt^nn — {gTrNNY^^ — 13.5 
during the numerical solving procedure of corresponding equations of motion. The typical 
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Figure 2.1: The radial dependence of chiral angle 9{r) - solid line and shape function 
p(r) = exp(— (T(r)) - the dotted lines for the soliton corresponding to the set of 
parameters U = SSMeF, e = 4.69 , Cg = {122.3MeV)^. 



9{r) and the shape function p{r) = exp{—a{r)) respectively. It is well known [27] that the 
Lagrangian (2.1) with Cg = (SOOMeV)'^ gives a large value for the soliton mass. Our precise 
analysis [29] show that the gluon condensate should be reduced up to the Cg = (122.3Me\^)^ 
to reproduce g^^N and the masses (Ma - = 293MeV, = 95A.5MeV) [29]. 

In this case the calculated value of axial coupling constant qa, which is given by the 
following expression: 

-TT rco 

g^^- — j^ dxx''{p\e' + S2lx)+A:[s2{e'^ + d)lx + 2e'<]\} (2.18) 

is equal to 1.01 i.e. much better then that in OSM [12] {qa^^^ = 0.61). It is clear from 
Eq. (2.18) that qa is very sensitive to the Skyrme parameter e which is connected to Qt^nn 
through the equations (2.6) - (2.8), since Goldberger Treiman relation does not work due to 

Thus by using the following optimal set of parameters: /^r = 93MeV , e — 4.69 , Cg ~ 
{122.3MeV)^ in Eqs. (2. 16), (2. 17) we calculated the scalar - isoscalar part of central nuclcon 
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Figure 2.2: Central isospin- independent potential as a function of internucleon dis- 
tance r calculated from Eqs. (2. 16), (2. 17) . The dashed curve is the corresponding 
interaction component in the "Paris" potential [41]. 

interaction component of the realistic phenomenological "Paris" potential [41], displayed 
there with the dashed line. It is seen from the figure that the attraction is shifted to larger 
separations and much stronger. The similar result was obtained in ref. [42] within the 
framework of realistic pseudoscalar - vector chiral Lagrangian where the a - meson exchange 
was simulated by a correlated two-pion exchange between the two solitons. 

To discuss this effect in a more detail we bring in figures Fig.2.3 (a) and Fig.2.3(b) the net 
contributions from each term. It is well known that in the OSM the second order derivative 
term £2 does not contribute to the central potential at all [p{x) = 1 in Eqs. (2.16), (2.17)], 
whereas the contribution from £4^^ is strongly repulsive (see Fig.2.3(b)). After the inclusion 
of scalar dilaton this repulsion may be compensated by an attractive contribution arising 
from the kinetic term £2 [see Fig.2.3(a)] in the region r > 1.5/m. As to the contributions 
from ^T^g]^ and W^j terms they are not large, especially, at short and medium distances [see 
Fig.2.3(a) and Fig.2.3(b)]. 

Note that the relevance of product approximation in the short range region of the inter- 
action may be doubtful. However, the considering attraction appears at separations greater 
than 1.2/m where the approximation is valid. 




Figure 2.3: The net contributions to the central NN potential given by the Eqs. 
(2.16), (2.17). The V2, Vgb, V4, and Vw curves are the contributions from -^^i^sb' 
C^Q^ and terms respectively. 

We conclude that the model we have presented here leads to a correct values of f^^, Qt^nn, 
MjvA and gives a desired attraction in the central part of NN interactions. On the other 
side, by comparing this result with a result of Yabu et al. [25] where no attempt had been 
made to optimize the parameters we may affirm that the attraction appears in a natural 
way in both cases regardless the origin of dilaton. One may expect that the inclusion of 
cu— meson coupling term to the Lagrangian [43] and taking into account of deformation 
effects [44, 45, 46] will give a better description of NN and NN interactions.That will be 
the subject of the next chapter. 



Chapter 3 



The Deformation of the Interacting Nucleon in the 

Skyrme Model. ^ 

3.1 Introduction 

The possible changes in the radius of a nucleon in interaction with another nucleon was 
investigated in the skyrmion model by Kalbermann et. al. [49]. It was found that there is a 
swelling of the nucleon at intermediate distances between the nucleons which were assumed to 
preserve their spherical shape. On the other hand Hajduk and Schwesinger [50] considering 
the skyrmion to be soft for the deformation showed that there are several deformed states 
of rotating skyrmions. In particular, besides the ground state of a spherical shape there 
exist a degenerate doublet exotic states with the same quantum numbers of the nucleon 
{s = t = 1/2) of oblate and prolate shapes. One may therefore wonder whether the nucleon 
can change its shape under the action of strong interactions. 

Certainly, several studies have been made of the deformation effects on the Skyrmions 
[44, 45, 46], but these were carried out to obtain a better understanding of the possible 
sources of attraction in the central nucleon - nucleon potential. As a result it was shown 
that the deformation effect is very important and may reduce the central repulsion by about 
40%. 

In the previous chapter the problem of missing central attraction in the A^A^ - interaction 
has been investigated within the model of Andrianov and Novozhilov [27] which starts with 
the Skyrme Lagrangian supplimented by a dilaton scalar field a{r) so as to satisfy the QCD 
trace anomaly constraint. It has been concluded that this model gives the desired attraction 
in the central part of the A^A^ interaction. The resulting central scalar - isoscalar part of the 
potential is in qualitative agreement with the phenomenological one. 

In the present chapter we shall concentrate on the effects of modification of the shape 
of a nucleon at a quantitative level using skyrmions. We shall use the same Lagrangian 
(2.1)- (2.4) [27, 51] which has been used in the previous chapter. 

3.2 The interaction of deformed nucleons 

As it was shown in the previous chapter, the Lagrangian CAN{U,a) is a generalization 
of the well - known original Skyrme model [12, 13] and takes into account the conformal 
anomaly of the QCD. 

We have also shown that the Lagrangian produces in a natural and transparent way the 
intermediate - range attraction in the central part of baryon - baryon interaction even in the 

^The present chapter is based on following articles by the author and his collaborators: [47, 48] 



product approximation [30]. Here we let the interacting nucleons deform. Assume that both 
skyrmions deform into an eUipsoidal shape as they come close together. To describe this we 
write the chiral field U and the dilaton field a as the nonspherical hedgehog form given by : 



Uo{r) = exp {ifq9{q)) ,^ 
a{r) = a{q) 

where the spatial vector g'has the components d^x^ d^yfizZ and q is the unit vector: q = q/q 
with the deformation parameters 6^ and 5z- The profile functions 9 and a are assumed to 
be the solutions of the Euler - Lagrange equations in the spherical case given in (2.6). The 
ansatz, in Eq. (3.1), leads to a modification of the static mass, M^, and the moment of 
inertia, A^, of the Skyrmion 

M*H = [AM2 + BM^, + M^,6 + Mw]h , . 

where r] = 616^, A = {25l + 5^)/3, B = 6l{5l + 25l)/3 and Mi and denote the relevant 
contributions from £i term in Eq.(2.1) for the spherical case 6x = = 1. As we are mainly 
interested in the region where the medium range attraction takes place - ~ 1.25/m ( 
typical separation between nucleons in nuclei) we restrict ourselves to the familiar product 

U = A.UoiX - ql2)AXMU^{X + ql2)A\ = U1U2 ^3 3^ 

P = P{X - q/2)p{X + q/2) = pips , P = cxp(-cT) 

where Ai, A2 are the collective coordinates of skyrmions to describe their rotational motion, 
q is the vector along z axis: qx = ^, Qy = ^, Qz = Q = RSz and R is the distance between 
skyrmions. The static skyrmion - skyrmion potential is defined by 

V(R, 5x, Sz)^-J dX[jC(UiU2, P1P2) - jC(Ui, pi) - jC(U2, P2)] (3.4) 

The application of the usual projection methods developed in [44, 45, 46, 52] to Eq.s 
(3.1) — (3.4) yields the following representation for the central scalar - isoscalar part of the 
nucleon - nucleon interaction : 

Vc{R, Sx, Sz) = -[Vxsbiq) + Vw{q) + Sl{V2{q) + SlVM)+ ^3 
+(Sl-5l){V2'^h<l) + SlV,t'm 

where the terms V'^'^-^ and V^^^ are the net contributions from the deformation effect of the 
terms £2 £4^^ in Eq.s (2.1)-(2.3) respectively. They are given by 

e Jo 

X /' dt[pi\p2' - l)(/3i + Y^tl) + 4>ltl{p2^ - 1)+ 

+P1P20102 + (1 ^ 2)1, 



3e JO 

X dt[{Y^xl + 3/?i)(r2i2 + 1^2) - Y1Y2ZZ, 
-/3i/32-2/32yii? + (1^2)]. 



with 

p{x) — exp {—a{x)), (f){x) — dp/dx, 

(5{x) = (sin e{x) jxf, Y{x) = \d'''{x) - /^(x)]/^^, 

Z = - gV4, Z3 = xH-" - q^/A, q = d^ref^, (3.7) 

^2 ^ = g2/4 _^ 3,2^2 _^ ^^^^ ^2^^ ^ ^2^^ + a;2(^;L _ ^2^^ 
Pi = p{Xi) (pi = (f){Xi) Yi = Y{Xi)J% = (3{Xi). 

Other terms in Eq. (3.5) don't include deformation effects exphcitly. They are brought 
in previous chapter. We shall calculate the deformation parameters 5a; (r) and 5z{r) by 

220 n 1 




0.50 0.7b 1.00 1.2b l.bb 1 . /b 2.00 2.2b 2.b0 

R( f m) 

Figure 3.4: Central isospin - independent potential in Eq.(3.5) calculated for the 
cases with the dilaton field (solid curve) and without one (dashed curve) as a function 
of the internucleon distance R . The dotted curve is the corresponding interaction 
component in the "Paris" potential [41]. 

minimizing the total static energy of two nucleon system at each separation R using Eq.s 
(3.1) - (3.5). 

The resulting values of the parameters will be used to study the changes in the shape 
of the nucleon. Now we illustrate this procedure for the case of the isoscalar mean square 
radius and the annronriate intrinsic nuadrunole moment. The normalized isoscalar mean 



square radius along each axis may be defined by 



^ . _ JdrrjBo(r) 
^^^^^=0" JdfBoif) ^^-^^ 

where i — x,y,z and So(r) is the baryon charge distribution 

Mf1-^e^''T^[LiLjL,]. (3.9) 

The inclusion of the deformation in a simple way: ri^qi/6i as in Eq. (3.1) yields the 
following relation between the radius of a free spherical nucleon (r^)j^Q and a deformed 

one: {r.i^)*j^Q = {i"'^) i=o- Therefore the appropriate quadrupole moment characterizing the 

shape of the distribution of baryon matter is compared to that of an ellipsoid with axis 1/Sz 
and l/Sx'. 

Qi=o = 5{r,ri=o - {rri=o = 2(r^),=o( V'^' - V^') • (3-10) 

The explicit formulas for Qj=i defined by Qi=i = 3{rz'^)*j^^ — {r^)*j^i are rather compli- 
cated and may by found elsewhere [48] . 



3.3 Results and discussions 

In the numerical calculations we consider the following two cases: the Lagrangian with the 
dilaton and the pure Skyrme model when o" = in Eq. (2.1). In both cases the parameters 

e and m^r were fixed at the values: = 93MeV, e = 27r, = 139MeV. For the gluon 
condensate we use Cg = (283MeV)^ as obtained from lattice QCD calculations [53]. The 
mass of the scalar meson, rria, defined by = ^2Cg/2/^ is then GlOMeV. 

This set of parameters produce the following static properties of the nucleon: Mjv = 
1054Mey, QA = 0.65, {r^)y^Q = 0.38/m and {r^^/^^ = 0.66/m in the dilaton case. No 
attempt is made here to search for a realistic set of parameters since our interest was mainly 
in the link between properties of A^A^ interaction and the shape of the nucleon. 

In Fig. 3.4 the central scalar - isoscalar part of the nucleon - nucleon interaction has 
been presented for both cases including deformation effects. For a comparison of the corre- 
sponding interaction component with the realistic phenomenological "Paris" potential [41] 
is also displayed here ( the dotted line in Fig. 3.4). It is clear that the Lagrangian with the 
dilaton field is able to describe the nucleon - nucleon interaction in the intermediate region 
quite well. 

We now turn to the changes in the shapes of interacting nucleons. The intrinsic quadrupole 
moments Qi=o and Qi=i are shown in Fig.3.5 and Fig. 3. 6 respectively. In the pure Skyrme 
model when there is no attraction [V2 = V^*^^ = in Eq.(3.5)] between skyrmions it becomes 
oblate (dashed lines in the figures 3.5, 3.6) due to the strong repulsion caused by the 
terms in Eq. (3.5). The inclusion of the dilaton leads to the following qualitative picture: 
At large separations skyrmion is obviously in a spherical shape, becomes prolate at the in- 
termediate region and deforms to an oblate shape at small distances where the repulsion 
dominates. 

As the nucleons approach each other they change shapes from prolate into oblate at 
R ~ 1.2/m. Comparing figures Fig.3.5 and Fig. 3. 6 it may be noticed that the isoscalar in- 
trinsic quadrupole moment (5/=o is much smaller than the isovector one Qi=\ at intermediate 
separations. The intrinsic quadrupole moment of the proton defined by: 
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Figure 3.5: The radial dependence of the isoscalar quadrupole moment Qi=o of the 
nucleon.The sohd and dashed hnes are obtained in the case with dilaton and pure 
Skyrme model respectively. 

reaches a maximum value of Qp = 0.016/m^ at r ~ 1.5/m. The corresponding excentricity 
e (see Fig. 3.7) reaches the maximum value e ~ 0.02. Note that for deutron e^i ~ 0.12, due 
to the small mixing of D - state. In this sense, one may conclude that the shape of a nucleon 
in the nuclei is "more spherical" than the deutron. We expect new data from high - energy 
electron scattering on nuclei to make the situation clear. 

As a concluding remark we have to underline that the deformed states of oblate (prolate) 
shapes may not necessarily belong to the K = 1 band found in rcf. [50] since for a strongly 
deformed system the quantization procedure used here needs some modifications. It is a 
generalization of the well - known original Skyrme model [12, 13] and takes into account the 
conformal anomaly of QCD. 



0. 04 




- . 24 \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ 

0.5 1.0 1.5 2.0 2.5 

R( f m) 



Figure 3.6: The same as in Fig. 3.5 but for the isovector quadrupole moment Qi=i. 




Figure 3.7: The excentricity of the elhpsoidally deformed nucleon defined as £ = 
(a^ — 6^) / (a^ + 6^) , where a and b are the tranverse and conjugate axes of the eUipse. 



Chapter 4 



The nucleon— nucleon interaction and properties of the 
nucleon in a -Kpuj soliton model including a dilaton field 

with anomalous dimension ^ 

4.1 Introduction 

Recently Furnstahl, Tang and Serot (FTS) [55] have proposed a new model for nuclear 
matter and finite nuclei that realizes QCD symmetries such as chiral symmetry, broken scale 
invariance and the phenomenology of vector meson dominance. An important feature of this 
approach is the inclusion of light scalar degrees of freedom, which are given an anomalous 
scale dimension. The vacuum dynamics of QCD is constrained by the trace anomaly and 
related low-energy theorems of QCD. The scalar-isoscalar sector of the theory is divided into 
a low mass part that is adequately described by a scalar meson (quarkonium) with anomalous 
dimension and a high mass part (gluonium), that can be "integrated out", leading to various 
couplings among the remaining fields. The application of the model to the properties of 
nuclear matter as well as finite nuclei gave a satisfactory description. Further developments 
of the model [56, 57] showed that the light scalar related to the trace anomaly can play a 
significant role not only in the description of bound nucleons but also in the description of 
heavy-ion collisions. It was also shown that the anomalous cannot be due to an effect of 
nuclear density on the trace anomaly of QCD [56] . 

Here a natural question arises: What is the role of this light quarkonium in the description 
of the properties of a single nucleon, when it is taken into account in topological nonlinear 
chiral soliton models, which are similar to the FTS effective Lagrangian on the single nucleon 
level? In the present Chapter we introduce a dilaton field with an anomalous dimension into 
the 7rpa;-model proposed in ref. [58] and investigate some properties of single nucleon which 
emerges as a soliton in the sector with baryon number one {B — 1). Note that, analysis of 
pion - nucleon scattering phase shifts made in refs. [59, 60] indicate the existence of the light 
scalar meson. 

It is well known that a scalar-isoscalar meson, the sigma, plays an important role in the 
nucleon-nucleon (NN) interaction especially within one-boson-exchange (OBE) models [61, 
62]. Note also that, the missing medium range attraction was a long standing puzzle in 
Skyrmc like models. In Chapter 2 we have shown that explicit inclusion of a scalar meson 
into the Skyrme model produces in a natural way the desired attraction. However, OBE 
model includes vector mesons also. Therefore, it would be quite interesting to investigate 
the central part of the NN interaction when the light scalar (a), p and a;-mesons are taken 
into account explicitly. This is what will be considered in present Chapter. In particular, it 
is important that when one is to properly describe the intermediate range attraction in the 

^The present chapter is based on following articles by the author and his collaborators: [29, 54] 



central NN interaction, the successful description of the single nucleon properties within the 
Trpo; model should not be destroyed. 

We also note that in an soliton approach with explicit regulated two-pion loop graphs 
one is able to get the proper intermediate range attraction. In that case, however, one does 
not stay within a simple OBE approach any more (as done here) and also needs to calculate 
the modifications of the isovector two-pion exchange to the p and so on. For comparison, we 
mention that recent developments of the original Bonn OBE potential performed at Jiilich 
also include multi-meson exchanges leading to a renormalization of various interactions, 
couplings and cut-off parameters [63, 64]. The model we investigate is related closely to the 
OBE approximation of the NN force. 



4.2 The npujG model 

Including a a-meson by means of the scale invariance and trace anomaly of QCD into 
the 7rpa;-model [65] can be done in terms of the following chiral Lagrangian of the coupled 
npuju system, 

C = ^ d^ad^'a - Tr L,L^^ - Tr [I, + r, + igfp^+ 

3 1 f^m^e-^"/'^ 
+igu;,r + i^gu^.B^ - -^i^,.^'" + m'") + ^-^^ TV {U - 1)- (4 ^) 

_^J«[l_e--/^(l^ + l)], 



16 ^ ^ d 

where the pion fields are parametrized in terms of [/ = exp {ir ■ tt / f^r) and ^ = \/U, 
left /right-handed currents are given by = U^d^U, = i~^dfj,^, r^j, = Cd^^^, and the 
pertinent vector meson (p, a;) field strength tensors are p^i, = d^p^ — d^Pj^ + g[Pn X Pv] 
and cu/j^i, — d/j^cui, — dvOJ^. Furthermore, the topological baryon number current is given by 

In Eq.(4.1), 5*0 is the vacuum expectation value of scalar field in free space matter, 
is the pion decay constant (/^r = 93 MeV) and g = gp^^^ is determined through the 
KSFR relation g = m/\/2/7r. The model assumes the masses of p and cu mesons to be 
equal, nip = = m. The mass of the a is related to the gluon condensate in the usual 
way [55, 66, 67] rua — 2yJ^g/{dSo), where d is the scale dimension of scalar field {d > 1). 
Being "mapped" onto the states of a nucleon, the Lagrangian Eq.(4.1) will be similar to the 
FTS effective Lagrangian. 

Nuclcons arise as soliton solutions from the Lagrangian Eq.(4.1) in the sector with baryon 
number i? = 1 as it has been explained in Chapter I. To construct them one goes through 
a two step procedure. First, one finds the classical sohton which has neither good spin nor 
good isospin. Then an adiabatic rotation of the soliton is performed and it is quantized 
collectively.^ The classical soliton follows from Eq.(4.1) by virtue of a spherical symmetrical 
ansatze for the meson fields: 

G(r) 

U{r) = exp(iff6l(r)), pf = Siakrk—^, ^u,{r) = oj{r)5^o, cr{r) = a{r) . (4.2) 

In what follows we call 0{r), G{r), 00 {r), and o"(r) the pion-, p-, uj-, and o"-meson profile 
functions, respectively. The pertinent boundary conditions to ensure baryon number one and 



finite energy are, ^(0) = 7r,G(0) = -2,uj'{0) = a'{0) = 0,^(oo) = G'(oo) = a;(oo) = ct(oo) = 
0. To project out baryonic states of good spin and isospin, we perform a time-independent 
SU(2) rotation 

U{f,t) = A{t)U{f)A+{t), e(r,t) = Am{f)A+{t) 

a(f, t)^a(r), uj (f , t)^^ [Kf] 

2 (4-3) 

T-p^{r,t) = -A{t)f- {KUr)+rk-ri2{r))A+{t), 

f ■ p^{r,t) = A{t)f ■ p^{r)A+{t) 

with 2K the angular frequency of the spinning mode of sohton, ir ■ K — A'^A. This leads 
to the time-dependent Lagrange function 

C{t) = j dfC^ -Mh{9, G, u, a) + A{9, G, u, a, 0, ^i, 6)Tr (ii+) . (4.4) 

Minimizing the classical mass MH{0,G,uj,a) leads to the coupled differential equations for 
6, G, uo and a subject to the aforementioned boundary conditions. In the spirit of the large 
A^c^expansion, one then extremizes the moment of inertia A(^, G, a;, cr, 0, .^i, (^2) which gives 
the coupled differential equations for ^1 , ^2 and in the presence of the background profiles 
9, G, (jj and a. The pertinent boundary conditions are (/)(0) = 0(oo) = 0, ^i(O) = Ci(oo) — 0; 
^^(0) = ^2(00) = 0, 2{i(0) +{2(0) = 2. The masses of nuclcon Mjv and the mass of A , Ma, 
are then given by Mn = Mh + 3/8 A and Ma = + 3/ 15 A. 

The electromagnetic form factors obtained in the usual way [65] are: 

G|(g 2) = J MqrMr)e-''^/'r'dr, 

^ 



(4.5) 



GUr) - ^ 7jo(?r) ||[44 + (1 + 2c)6 + 6]e-^^/'^ + '-^Vdr, 
^ 

where s — sin (^), c = cos(^) and S2 = sin(^/2). The normalization is ^1(0) = G^(0) = 1/2. 
Similarly, meson-nucleon vertex form factors may be calculated [68, 69]. In the Breit frame 



one has: 
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where the "electric" and "magnetic" vector meson-nucleon form factors are connected to 
the Dirac Fi{t) and Pauh form factors -F2(^) through the following relations: 



GUt) = Fl{t)+tFm/^Mj,, 
G^^{t)=Fl{t) + Fm, (^ = P,a;). 



(4.7) 



4.3 Results and discussions. 



4.3. 1 Static and electromagnetic properties of the nucleon. 

Using the formulas given above we have calculated static and electromagnetic properties 
of nucleon. As can be seen from Eq.(4.1), the Lagrangian has no free parameters in the Trpcu 
sector. So, in actual calculations the parameters mTr, m, f^, are fixed at their emperical values, 
= 138 MeV, m = rrip = = 770 MeV, = 93 MeV, g = m/V2U = 5.85. In the 
(7-meson sector there are in general three free parameters: m^r, Sq and the scale dimension 
d. The latter has been well studied for nuclear matter calculations [55, 56]. In particular, 
it was shown that for d>2 the much debated Brown-Rho (BR) scaling may be recovered. 
Therefore, assuming that there is no dependence of d on the density, wc shall use the best 
value d = 2.6 found in refs. [55, 56]. The values for 5*0 - were found to be 5*0 = 90.6^95.6 MeV 
[55] . So we put So = fn = 93 MeV. The mass of the a or equivalently the gluon condensate 
Cg = m?^d'^Sl/A = m?^d'^fl/A is uncertain. We thus consider two cases: rricr — 550 MeV 
and Too- — 720 MeV in accordance with recent tttt phase shift analyses [59, 60] and with 
OBE values. We stress again that the precise nature of such a scalar-isoscalar field is not 
relevant here, only that it should not be a pure gluonium state. A summary of static nucleon 
properties obtained in both cases, i.e. with m„ = 550 MeV and to^. = 720 MeV, is given in 
Table 4.2. 

One immediately observes that the nucleon mass is again overestimated. This may not 
be regarded as a deficiency, since it is known that quantum fluctuations tend to decrease the 
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Table 4.2: Baryon properties in the irpu and irpua models 





TVpUJ 


Tvpooa 


Tvpooa 


Exp. 


[MeV] 


— 


550 


720 




[MeV] 


— 


258 


295 


300^400 


[MeV] 





119 


121 





Mn [MeV] 


1560 


1492 


1511 


939 


A [fm] 


0.88 


0.88 


0.88 


— 


Ma - Mjv [MeV] 


344 


350 


350 


293 


rn = {rlY'^ [fm] 


0.5 


0.5 


0.5 


-0.5 


(r|)y2 [fm] 


0.92 


0.94 


0.94 


0.86± 0.01 


(rl)n [fm^] 


-0.20 


-0.16 


-0.16 


-0.119±0.004 


(rl.)y'[fm] 


0.84 


0.85 


0.85 


0.86± 0.06 


(rl.}^/^ [fai] 


0.85 


0.85 


0.85 


0.88± 0.07 


Hp [n.m.] 


3.34 


3.33 


3.33 


2.79 


/i„ [n.m.] 


-2.58 


-2.53 


-2.53 


-1.91 


I)"p/)"n| 


1.29 


1.30 


1.30 


1.46 


9A 


0.88 


0.95 


0.95 


1.26±0.006 




0.63 


0.66 


0.66 


0.65± 0.07 



To estimate the influence of the cr-meson we also show the results given by minimal 
version of Trpcu model. As can be seen from Table 4.2, the inclusion of a light sigma meson 
into the basic npuj model just slightly changes the nucleon mass and its electromagnetic 
properties. This may be explained by the fact that the role of intermediate scalar-isoscalar 
meson in gamma-nucleon interactions is negligible. In contrast, the presence of the sigma- 
meson leads to an enhancement of axial coupling constant as it was first observed in ref . [27] . 
Although, the physical mechanism of this change is not clear, it may be understood as mainly 
due to modification of meson profiles. This is in marked contrast to the inclusion of pion 
loop effects, which tend to lower the axial coupling even further [70]. In the present model 

= 0.88 for the Tipuj and Qa = 0.95 for the npuja model, respectively. One may expect that 
an appropriate inclusion of a meson into a more complete version of the basic Trpa; model 
might give the desired value = 1-26. 

4.3. 2 Meson— nucleon form factors. 

One of the usual ways to calculate the meson-nucleon interaction potential within topo- 
logical soliton models is the so-called product ansatz [25]. Within this approximation, 
the two-Skyrmion potential as a function of the relative angles of orientation between the 
Skyrmions has a compact form, and the extraction of the NN potential by projection onto 
asymptotic two-nucleon states is straightforward. This procedure gives only three nonva- 
nishing channels: the central, spin-spin and tensor potential. At large and intermediate 
distances, the latter two compare well with e.g. the phenomenological Paris potential [41] . 
The major inadequacy found in such type of calculations is the lack of an intermediate range 
attraction in central potential. Although many remedies have been proposed, this result 
may not be genuine for Skyrme like models. In fact, the product ansatz, which is not a 
solution to the equations of motion, can only be considered accurate at large distances, and 
the failure of these calculations to reproduce the central range attraction may simply be 



the failure of the product approximation to provide an adequate approximation to the exact 
solution.^ Although the lack of central attraction may be recovered by the inclusion of a 
scalar-isoscalar meson, the inherent ansatz dependence of the trial configuration remains as 
a major shortcoming of product approximation [52]. 

On the other hand there is another natural way which was first used by Holzwarth 
and Machleidt [72]. They proposed to calculate Vjvtv within OBE model taking coupling 
constants and meson-nucleon form factors from a microscopical model such as the Cloudy 
Bag model or the Skyrme model. It was shown that the Skyrme form factor is a soft pion 
form factor that is compatible with the vrA^ and NN systems. We shall use this strategy to 
investigate the NN potential within present model. 

The meson nucleon form factors given by well known procedure, proposed first by Co- 
hen [69] are given in Eq.(4.6). Although they were derived in a microscopical and consequent 
way, these form factors could not be directly used in standard OBE schemes. The reason is 
that the OBE schemes [62] in momentum space use form factors defined for fields propagat- 
ing on a fiat metric, whereas the definition of form factors in Eq.(4.6) involve a nontrivial 
metric. Hence, before using the latters in OBE scheme one should modify the procedure 
by redifining meson fields. The modification for pion nucleon form factors in Trpcu model is 
clearly outlined in refs. [72, 73]. Now, applying this procedure in the Lagrangian in Eq. 
(4.1) we get the following nNN form factor: 

G-^-q^) = ^l^MEA^^2^rnl) Jji(gr)0wJ(O sin(^)r2dr = 



— / — b — -\—2F [r] — rF (r) -I — + rmz.F{r)\r dr, 

3 J qr r 
^ 



(4.8) 



where 

Mt = [1 + 2tan2(^/2)] e-^'^/^ 

/ (4-9) 

F{r) = ^3 + cos2(^(r)) - 4 cos(^(r))e-'^/<^. 

The infiuence of metric factor Mt to the pion-nucleon form factor is illustrated in Fig. 4.8. 

It is seen that, without the inclusion of Mt the form factor is softer than in OBE models 
(the dashed fine in Fig. 4.8), while its inclusion via Eqs. (4.8), (4.9) gives a behavior closer 
to OBE models. In fact, a monopole approximation at small q^ — t oi the normalised form 
factor G^(t)/G^(0) ^ A2/(A2 - t) gives A^ = 860 MeV and A^ = 1100 MeV for Mr = 1 
and Mt 7^ 1 respectively, compared to its emperical fit : A^^^ = 1300 MeV (dotted line in 
Fig. 4.8). Note, however, that our result for A^r including the metric is in line with recent 
coupled-channel calculations of the Jiilich group [74]. There, a monopole form factor with 
A^ ~ 800 MeV is obtained. 

Introducing a fiat metric requires a canonical form for the kinetic part of the Lagrangian, 
which determines the dynamics of the field fiuctation. The kinetic term of the scalar meson 
in Eq. (4.1) 

C^''^ = S^e~^''/%ad^a/2 (4.10) 
can be easily rewriten in a usual way: 

= d^ad''a/2 (4.11) 



^An early study giving credit to this line of reasoning can be found in ref.[71]. 




Figure 4.8: The normalized ttNN form factor in the npuoa model {nia- = 720 MeV). 
The solid line represents the form factor when the metric factor is included (Eq. 
(4.8)), while the dashed hne gives the result with no metric factor as in Eq. (4.13). 
The dotted hne is a monopole form factor with A^r = 1300 MeV. 



by the following redifinition of the basic sigma field: 

a(r) = Sod[l - g-'^^'')/'^] 



(4.12) 



Now the new field a may be identified with the real sigma field. Clearly this redifinition 
docs not change the nuclcons static properties given in Table 4.2. Note also that, using the 
above redinition in the last term of Eq. (4.1), one may easily conclude that = ni^r- The 
appropiate sigma-nucleon form factor is given by 
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mla 



r^dr, 



(4.13) 



and may be used in OBE models. We have not introduced any metric factors in the form 
factors of the heavier mesons since these should play a lesser role than in the case of the 
pion. 



For small values of the squared four-momentum transfer t each form factor can be 
parametrized in monopole form: 

G,{t)=g,{A^-m^)/{A^-m^) (4.14) 

{i — TT, p, a;,(7). We present in Table 4.3 the range parameters (cut-offs) and the couphng 
constants of the resulting meson-nucleon dynamics. 

Table 4.3: Meson-nucleon coupling constants and cut-off parameters of meson- 
nucleon form factors. The Ai{i = tt, p, ua) are cutoff parameters in equivalent 
monopole fits 1/(1 — t/A|) to the normalized form factors GiNN{t)/GiNN{0) 
around t = 0. The empirical values are from OBE potential fit [62] 





npu 


7rpc<j(T(550) 


7rpLua{720) 


OBE/Emp. 




14.74 


13.97 


14.17 


13.53 


GaNN{Q) 




6.2 


6.19 


9.1 (12.41) 


Ff(0) 

Fm 

Fm/Fm 


2.55 


2.76 


2.68 


2.24 


14.33 


15.01 


14.67 


13.7 


5.6 


5.43 


5.47 


6.1 


Frio) 


8.78 


10.73 


10.15 


11.7 


F-(0) 


-2.15 


-2.78 


-2.65 





F-(0)/Fr(0) 


-0.24 


-0.25 


-0.26 





A. (GeV) 


1.2 


1.1 


1.1 


1.3^2.0 


A. (GeV) 
A1 (GeV) 
A^ (GeV) 




0.59 


0.60 


1.3^2.0 


0.62 


0.63 


0.63 


1.3 


0.92 


0.92 


0.92 


1.3 


A^' (GeV) 


0.95 


0.89 


0.91 


1.5 


A- (GeV) 


1.12 


0.86 


0.89 





One can see that there is no interference between cr-meson nucleon and e. g. the pion- 
nucleon coupling constants. In other words, the inclusion of cr-meson does not significantly 
affect meson-nucleon form factors that had been given by the Trpcu model. As it is seen 
from Table 4.3 the values for meson-nucleon coupling constants are close to their emperical 
values (in some cases obtained by OBE model fits). This is one of the main advantages of 
the inclusion of a scalar-isoscalar meson as done in the present approach. 

The corollary of the present model is that it gives significant information on the a- 
nucleon interaction. As it is seen from the Table 4.3, the value for QaNN and the cut-off 
parameter of sigma -nucleon vertex Aq- are smaller than their OBE prediction K^^^ ^ 
1300 ^ 2000 MeV. This contrast is evidently seen from the Fig. 4.9, where G"" (t) / G"" (0) for 
two cases : 171^ = 720MeV and = 550MeV is presented with the solid and dashed lines 
respectively. The band enclosed by the dotted lines refers to the OBE monople form factor 
with A^^^ = 1300 . . . 2000 MeV. One can conclude that the present model gives a softer 
aNN form factor than it obtained by OBE. As it had been noticed before, the t-plane for 
each form factor has a cut along the positive real axis extending from t = to to (yo. The cut 
for the cr-nucleon vertex function starts at to = 4m^^ reflecting the kinematical threshold 
for the (7 — > TTTT channel. More precisely this result follows from the asymptotic behavior of 
the meson profiles: For r — > oo, we have 9{r) ~ exp(— m7rr)/m^r, a{r) ~ ^^(r), which are 
derived from the equations of motion. 
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Figure 4.9: The sigma - nucleoli form factor GaNN{(l'^)/GaNN{0)- The dashed and 
solid lines are for rua — 550 MeV and rua — 720 MeV, respectively. Typical OBE 
monopole fits with A = 1.3 ... 2 GeV are shown by the band enclosed by the dotted 
lines. 



4.3. 3 Nucleon - Nucleon central interaction 

Once the vertex function of the corresponding meson-nucleon interaction has been found, 

its appropriate contribution to the A^A^ interaction may be easily calculated by using well 
known techniques from OBE. The detailed formulas are given elsewhere [42, 62]. In partic- 
ular, the contribution of the cr-meson exchange to the central potential is given by 



The central A^A^ potential in the T = 0, 5" = 1 state (the deuteron state) is presented in 
Fig. 4.10 in comparison with Paris potential. Our prediction is in good agreement with the 
emperical one. Note that the desired attraction in the central Vnn has before been obtained 
in the Trpu! model by means of two-meson exchange [42]. 

In conclusion it should be noted that we do not intend to describe (cover) all A^A^ phase 
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Figure 4.10: The central potential in the S — 1, T — state for npuj and npuja 
models (dashed and solid lines, respectively). No contribution from two-meson 
exchange has been taken into account. The dotted line corresponds to the Paris 
potential [41]. 



are usually included in OBE picture, the full model should also take in account e.g. NAp 
couplings. In addition the 2it exchange and its strong mixing with cr meson exchange (see 
e.g. ref.[75]) should be considered. Another reason possible which limits the accuracy oi NN 
phase shifts in the present model is that the crNN coupling is not sensitive to the mass of 
sigma (see Table 4.3) as it is in the OBE phenomenology. In fact, even when the 27r exchange 
is disregarded , the pure OBE model has to consider two types of sigmas with nearly the 
same masses but with quite different coupling constants. So, we refrain from performing 
direct calculations of NN phase shifts in the present model. Instead, we point out that the 
meson-nucleon form factors found in the present model could be useful in a wider context of 
calculations of nucleon-nucleon observables (phase shifts, deuteron properties etc) and may 
give more information on meson-nucleon and nucleon-nucleon dynamics. 



4.4 Summary 



an explicit light scalar-isoscalar meson field, which plays a central role in nuclear physics, 
based on the chiral symmetry and broken scale invariance of QCD. We have shown that for 
the single nucleon properties the successfuU description of the electromagnetic observables of 
the Tipuj model is not modified and even the value for the axial-vector coupling is somewhat 
improved. In the two-nucleon sector, this extended npuja Lagrangian leads to the correct 
intermediate range attraction in the central potential and a soft aNN formfactor for both 
values of sigma meson mass rua — 550 MeV and iria — 720 MeV. 



Chapter 5 



Density Dependence of Nucleon Properties in Nuclear 

Matter. ' 

5.1 Introduction 

It is well established that in - medium nucleon - nucleon (N-N) cross sections manifest 
in heavy ion collisions need not to be the same as their free space values [76, 77]. The 
origin of these changes must be reflected in modiflcations of the N-N interaction which is 
predominantly one-boson exchange. The parameters in meson exchange models are meson 
nucleon coupling constants and their masses. To avoid divergencies in loop integrals the 
meson - nucleon vertices are modified by form - factors which effectively provide cut-off 
parameters. All these parameters would change in the nuclear medium. 

Even for the free N-N interaction the boson exchange models use a phenomenological 
ansatz for vertex form factors. Recently Holzwarth and Machleidt [72] have shown that 
among the QCD inspired pion - nucleon form factors the Skyrme form factor [69] is the most 
preferable: it can describe well both the pion - nucleon and nucleon - nucleon systems. 

The purpose of this Chapter is to investigate the role of the medium in modifying the 
properties like meson nucleon coupling constant, meson- nucleon form factors and the axial 
vector form factor and the coupling constant. The calculations provide a consistent approach 
to modification of these properties in a modified Skyrme model that may be valid in a nuclear 
medium. The cr-meson is introduced as a dilaton field to satisfy scale invariance. Even though 
it is well known that scale invariance is badly broken, this provides a way to introduce the 
(T-meson which is so essential for an N-N interaction that is consistent with the two-body 
scattering data and the bound deuteron. 

5.2 The modified Skyrme Lagrangian 

In ref.s [78, 80] a medium modified Skyrme Lagrangian was proposed and applied to 
study the static properties of nucleons embedded in the nuclear medium. It gave a good 
description of changes in nucleon mass and its size in the medium. Here we shall outline 
the basic features of the Langrangian and then extend it by including the dilated a meson 
field in order to satisfy the scale invariance. It is well-known that scale invariance is badly 
broken, this symmetry is retained here to investigate its consequences. 

^The present chapter is based on following articles by the author and his collaborators: [78, 79, 80, 81, 82] 



Our basic assumption in modifying of the Skyrme Lagrangian 



^sk = T^Tr {d,U){d^U+) + ;^Tr [U+d,U, U+d^?- 



(5.1) 

-^Tr([/+-l)(t/-l), 



where and e are the parameters of the model, was that, in the lowest order expansion of 
the chiral field U — exp{2ir7r / F^^) fa 1 + 2i{r7r)/F„ — 2n'^/F^ + . . . the appropriate medium 
modified Skyrme Lagrangian should give the well known equation for the pion field [85] : 

df^d^Ti + {ml + fl)7r = 0, (5.2) 

where m^r is the pion mass and 11 is the self energy operator. This may be achieved simply 
by including 11 in the pion mass term of the Skyrme Lagrangian: ^ 

^Xsb = -^Tr [([/+ -!)(! + fl/ml){U - 1)]. (5.3) 

In general, the operator fl = fl^ + Hp in Eq. (5.3) acts both on the center of mass 
coordinate, R, and on the internal coordinate, r, of the Skyrmion. Only homogenous nuclear 
matter is considered. The translational invariance of the medium and that of the basic 
Skyrme Lagrangian (5.1) makes the coordinate dependence of 11 simpler: IT = n(r — R) 
which is also relevant for a moving Skyrmion. Further we assume that the Skyrmion is 
placed at the center of the nucleus i.e. R = 0. 

The P - wave part of the pion self energy IIa is dominated by the P33 resonance. A 
simple model used in practical calculations is the delta - hole model which concentrates on 
the pion - nucleon - delta interaction ignoring the nuclcon particle - hole excitations. In 
momentum space the IIa is given by Ilp{uj, k) = IIa(^^, k) = —k^xi^^ +fi'oX(^) ^)) [85] 
where g'^ is the Migdal parameter which accounts for the short range correlations, also known 
as the Ericson - Ericson - Lorentz - Lorenz effect. The pion susceptibility x has nearly linear 
dependence on the nuclear density p: 



where f-^NA ~ '^f-nNN ( I^nn/^'^ ~ 0.08) is the coupling constant, loa ~ P/2Ma + Ma — Mjv 
and b {b ^ is the range of the vertex form factor that is chosen to have a Gaussian 

form. For the pions bound in the nuclear matter cu is small < a; < m^r, so in the lowest 
order IIa has the form: 

UA{u;,k)^-PxA-^^, (5.5) 



where 



XA = 47rcop/(l + 47rt/», 

art = Xa/(1 + ^ng'^cop), (5.6) 

Co ^ Sf^j^jQmlAnuA ■ 

In the static case (a; — > 0) it coincides with the Kisslinger optical potential [85] which was 
used in ref . [80] . 



^Here and below the asterisk indicates the medium modified operators in quantities. 



Then in coordinate space the following symmetrized Lagrangian may be obtained 
- (f ) (^) - §(1 - x.)T.(W-)(W)+ 



32e2 



where m*^ = 771^^(1 + ns(p)/m^2), the effective pion mass arises from the S - wave part 
of the self energy Ils{p). Since the operators in nuclear matter do not satisfy the Lorcntz 
invariance there are two different effective coupling constants in the kinetic terms of Eq. 
(5.7) ^. The third term with "mixed" derivative in the Lagrangian in Eq. (5.7) clearly 
vanishes in free space {art = if p = 0) and contributes mainly to the moment of inertia, 
I, of the skyrmion. We shall consider this contribution in Sect. 5 of the present Chapter by 
estimating the NA mass splitting. 



5.3 Inclusion of scalar meson 

Recently arguments have been given in favor of the empirical evidence for a scalar- 
isoscalar meson, a. For example the isosinglet resonance with a mass mg. = 553.3 ± 0.5 
MeV and width F = 242.6 ± 1.2MeV^ [59, 60] found in the recent tttt phase shift analyses 
is believed to have the properties of the a meson that is essential for the N-N potential in 
OBE models to explain the NN-scattering data. 

However the a meson as a chiral partner of the pion was originally excluded from the 
Skyrme type Lagrangians. The only way of including a meson in the Lagrangian is by means 
of a dilaton field appropriate to the scale transformations — > e^x^. This procedure has 
been outlined in Chapter 1. 

Now a skyrmion imbedded in nuclear matter is considered. When a dilaton field is 
introduced into the effective Lagrangian the dilaton potential V{a) would be modified by 
a a field generated by the medium itself. While several attempts have been made [88] , the 
correct nature of this modification is still poorly understood. On the other hand it is well 
known that the gluon condensate Cg as well as the quark condensate (qq) decrease in the 
medium due to the partial restoration of the chiral symmetry. The present study is restricted 
to considering the medium modification of the dilaton potential, Eq. (1.42), by taking into 
account mainly the change of Cg i.e. Cg ^ C*. 

Thus putting all these considerations together the following Lagrangian is proposed for 
homogeneous nuclear medium in the static case 

^Ik = ^^Tr + J-^TY [L„ L,f + ^X^Tr {U+U^-2)- 

(5.8) 



7712 r<* 
8 ^ ^ 24 



i-x'-^(i-x1 



where xi^) = ap = 1 — xa, = U^diU . Note that the Skyrme parameter e coupled 

to the fourth derivative term remains unchanged since this term is related to the exchange 
of a very heavy p - meson with mass m* = rUp —>■ 00 [34]. 



There may be two alternative approaches to applying this Lagrangian in nuclear physics. 
In Quantum Hadrodynamic models, the a field plays the role of an external field modifying 
the properties of a soliton [88] or a bag [89] which moves in the background generated by 
the medium. In the present model this approach would mean that a = cr^R) and U = U{r) 
[88] that makes the scale invariance doubtful. In contrast, the mean field approximation is 
not used for the a field. Instead, the cr-field is strongly coupled to the nonlinear pion fields 
so as to generate the sohton . We assume that a = a{r — R) and U — U{r — R) ior a, moving 
skyrmion. In ref. [80] a similar approximation (Eq. (5.8) with a — 0) was used to estimate 
the medium modified static properties of the nucleon and found a well known behavior of the 
nuclcon mass M^/Mjv < 1 and its size i?^/i?jv > 1- In the next section we shall investigate 
in detail the dynamical properties of the meson - nucleon system. 



5.4 Meson - nucleon form factors and Goldberger - 
Treiman relation 

The semiclassical procedure for calculating the meson - nuclcon vertex form - factors in 
a topological chiral effective Lagrangian [68, 69] is well-known. In fact, the results of more 
accurate methods [90, 91], based on the correct quantization of the fluctuating chiral fields 
nearly coincide with the original result that was given by Cohen et al. [68] . Using the ansatz 
U{f,t) = A{t)Uo{f- R{t))A+{t) and defining the pion field as 

= -^Tr [TM{r- R)A+] (5.9) 
the following expressions are obtained 



3 J qr 
^ 

for the pion nucleon form factor and 

00 

G:^M = 27rF,(g 2 + mf) J Mqr)a{rydr 



(5.10) 



00 

27rF, J jo(qr)S,(r)dr 




(5.11) 



for the sigma nucleon form factor respectively. The details and the explicit expressions for 
the source functions S.,^{r) and Sa{r) are given in the Appendix A. Here we note that, in 
the chiral soliton models formulas for meson nucleon form - factors are mainly determined 
by the quantization scheme rather than by details of the Lagrangian. The latter manifests 
itself through equations of motion whose sohitions 9{r) and a{r) with spherically symmetric 
ansatz are Uo = cxp{irn6{r)), n = f/r, a{r) = a{r) that should be used in Eq.s (5.10), and 
(5.11). The effective mass of the nucleon is given by 

M*N^M*H + ^^ (5.12) 
where is the mass of the classical hedgehog soliton and /* is the moment of inertia of 



The effective pion decay constant /* sliould be understood before considering tlie Gold- 
berger - Treiman (GT) relation. Medium renormalized pion - decay constant /* can be 
naturally defined by the PCAC relation: 

VAa{x)=f:ml\a{x) (5.13) 

The medium modified axial coupling constant g\ measures the spin - isospin correlations 
in a nucleon, embedded in a medium and is defined as the expectation value of the space 
component of the axial current A]^ in the nucleon state at zero momentum transfer [83, 84] : 

]im{N{P')\AUr)\N{P)) = | hmG^^lg 2)(iV|^|iV) exp(^gr) (5.14) 

where q — {P' — P) and G^(g ^) is the axial form factor of nucleon, G\{Gi) = g\. Here Ui is 
the component of the nucleon spin. Due to the semiclassical quantization prescription, Eq. 
(5.9), the matrix element of the pion field evaluated between nucleon states is given by: 

(Ar(P')ka(r - R)\N{P)) = -^^^^P(^g^ J dxe-''^ sm{e) {N\aa{Tx)\N) (5.15) 
Evaluation of the matrix elements between nucleon states for both sides of Eq. (5.13) yields 

g\ = / sin(^)r=^dr (5.16) 



that was originally derived in [92] for a free particle. By comparing this equation with the 
expression for pion nucleon coupling constant: 5'*^7v — '-^^Afiv(9^)lg=o given by Eq. (5.10) 
the following medium modified Goldberger - Treiman relation is realized 

9:NNf: = 9*AM;, . (5.17) 

Although this relation has been suggested earlier in ref.s [93, 101] it has not been proved 
yet. On the other hand g\ and the axial form factor GAiq^) may be calculated directly from 
the Lagrangian in Eq. (5.8) in terms of the Noether currents [83, 84]. This gives 



oo 



GA{q^) = -47r / [jo{qr)A^{r) + ^-^A2{r)ydr, 
J qr 





with 



A(r) = ^ 
^ ' 8r 



0' ( „o 8(i 



(5.18) 



A2(r) = -Ai(r) + -(FXe-2<^+ ^2 



and 



oo 

g*A^-\j {e-''^('-)F>, (e' + 7) + ^ [(^" + + 2^'^] } ^'^^ (5-19) 



where d = sin^(^)/r^, S2 = sin(2^), and ap is defined in Eq. (5.8). The renormalized pion 
decay constant /* is obtained by combining the results given in Eq.s (5.17) and (5.19). 



Table 5.4: Density dependence of hadrons properties. All quantities are in MeV. 

p/po ml ml {C;y/' T;^,, ^M°^ 6M'^^ 

0.0 139.00 550.1 260.70 1413 251.2 283.7 283.7 

0.5 144.90 513.8 251.06 1271 88.7 200.1 238.1 

1.0 149.06 493.8 246.12 1157 34.6 161.9 205.4 



5.5 Renormalization of hadron masses 

Before going to a quantitative analyses of the medium effects we fix the following set of 
parameters for free space: = 2fl = ISQMeV, m^ = 139MeV, Cy = (260Me\/)^. This 
gives a good description of the sigma meson properties m^^ = 550MeV, Fg. = 251.2Mel^, 
that may be compared with their experimental values obtained from the recent tttt phase shift 
analysis [59, 60]. The Skyrme parameter e has been adjusted to reproduce the pion nucleon 
coupling constant: g-j^NN = 13.5 for e = 4.05. The well-established fixed parameters of the 
P - wave pion self energy in Eq. (5.5) are used in the pion sector: g'^ = 0.6, cq — 0.13m~^ 
[80, 85]. 

Now, the medium dependence of the input parameters may be considered. The possible 
renormalization of the Skyrme parameter e cannot be studied in the present approach unless 
the p meson is included in the Lagrangian explicitly. So we take e* — e. 

We adopt the following parametrization of m* 



ml = m^yl + Ils{p)/ml = m.„^l - ATxhopr]/ml , (5.20) 

where 77 = l+m^^/ Mj^ and is an effective S - wave ir — N scattering length. It is anticipated 
[80] that the results will not be too much sensitive to the value of bo. 

The only input parameter in the scalar meson sector is C*. The medium renormalization 
of the gluon condensate C*, in contrast with the renormalization of the quark condensate 
(qq) [94] and meson masses [95], is poorly known. However in the present approach (5.8) C* 
may be determined by m* through the equation 



SF^ml^N 



Various approaches [89, 93, 101] show that m*, has a linear density dependence. The fol- 
lowing parametrization 

^ = 1-0.12-^ (5.22) 
mo- Po 

is adopted here. It is consistent with the one obtained in the Quark coupling model (QCM) 

of framework [89] . 

The results for the static properties of hadrons are presented in Table 5.4. The second 
(m*) and the third (m*) columns of the table should be considered as input data for they 
were taken from other models [85, 89]. The medium renormalized gluon condensate C* is 
calculated from Eq.s (5.21) and (5.22) with Nf = 2, e = 8Nf/{33-2Nf). The change in the 
gluon condensate is small ~ 5% at normal nuclear matter density, po = 0.5m^. The stiffness 
of the gluon condensate as a consequence of the lack of scale invariance of QCD has been 
shown by Cohen [96] who found that the fourth root of the condensate might be altered by 
no more than 4% [94]. 

The main contribution to the ami vertex, and hence, to the decay width of a meson at 



the first term of the Lagrangian in Eq. (5.8). The Table 5.4 shows that the width Fo-^tttt is 
decreased significantly in the medium. This stimulates an interest to observe the a mesons 
in nuclei by experiments proposed recently by Kunihiro et el. [97]. 

Now medium effects on the mass of nucleon Mj^ and A^A mass splitting 5Mjva = ~ 
M]^ will be considered. In general, it is almost impossible to reproduce simultaneously the 
experimental values of masses and coupling constants within the Skyrme model even for 
a free particle. Since dynamics is the main interest, the set of parameters, was chosen so 
as to reproduce the pion nucleon coupling constant: Qt^nn — 13.5. It is clear from Table 
5.4 that the free space value of the nucleon mass is slightly large = 1413Me\^, 
whereas 6Mna is reproduced rather weU 6Mna = 28AM eV (SM^^ = 293MeV). The 
effective mass of the nucleon in normal nuclear matter density is decreased by a factor 
of M^/Mn — 0.82 which is in a good agreement with the estimates based on QCD sum rules 
M^fiQCD) = 680 ± SOMeV i.e. M^/Mn = 0.72 ± 0.09 [98]. We underline that the mass of 
the nucleon should be treated as the mass of the baryon which emerges as a soliton in the 
sector with baryon number one (B = 1). 

The study of NA mass splitting gives a chance to estimate the contribution from the 
"mixed derivative" term - the third term on the r.h.s of Eq. (5.7): 



Actually, owing to the canonical quantization SMna is related to the moment of inertia 
/ by: 5Mna = — = 3/2 (Jq + J^), where is the net contribution from Crt (clearly 
Irt = in free space). The explicit expressions for Iq and Irt are given in the Appendix A. In 
Table 5.4 are shown SM^a that has been calculated with the inclusion of Irt (denoted here 
SM^a) ^-^d without the inclusion of IrtiSM^A.)- the nuclear medium the Crt term leads 
to an enhancement of the moment of inertia decreasing 5Mjva significantly. Even without 
the term Crt in the Lagrangian the shift of SMj^a from its free value 6M'^^/6Mi\fA = 0.75 is 
larger than that obtained by Meissner [93] in a medium modified chiral soliton model based 
on the Brown Rho (BR) scahng law {5M*^J5Mna = 0.87). 

5.6 Renormalization of coupling constants and form - 
factors 

The axial - vector exchange currents must be considered in order to investigate the 
medium effects on qa and on the axial form factor G^(g^). However it is known that the 
bulk of exchange current effects arise from the A-hole contributions. Including such A-h 
effects would imply that bulk of the exchange currents effects are included in the effective 
qa- In a heavy nucleus it is meaningful to take these axial exchange operators into account 
as corrections to the effective axial current operator of a single nucleon Aq, = —g\&fa- So 
g*A in Eq. (5.19) may be considered as an effective axial coupling constant modified by the 
medium polarization and screening effects since we are considering an effective one body 
problem of a nucleon embedded in the nuclear medium. 

The second column of Table 5.5 displays a well known quenching behavior of qa that is 
mainly caused by a factor ctp = 1 — xa < 1 in the first term of Eq. (5.8). Note that the 
same set of input parameters F^, e, Cq (but without dilaton field a — gives the desired 
ratio g*A/gA — 0.8 [80]. The present calculations show that the inclusion of the scalar meson, 
which induces an additional attraction, prevents larger quenching: g\lgA = 0.9. 

In nuclei, a nucleon polarizes the medium in its vicinity. This leads to a screening 




(5.23) 



Table 5.5: Coupling constants and cut-off parameters at finite density. (All values are nor- 
malized to their free space ones.) 

P/PO 9*a/9A glNNhnNN 9*^NN/9aNN f^/U Kl K K/^a ^a/^a 

0.5 0.96 0.91 0.88 0.98 0.70 0.90 0.93 

1.0 0.92 0.80 0.78 0.94 0.56 0.84 0.55 



the screening mechanism may be described as being due to virtual Ah excitations that have 
been taken into account by the self - energy Ha term in Eq. (5.5). At normal nuclear matter 
density the renormalization of QnNN amounts to a reduction of 25% of the coupling strength. 
This is sufficient to explain the quenching of the Gamov - Teller strength in heavy nuclei [99] . 
Furthermore this is consistent with a general argument based on Ward-Takahashi relations 
[100]. 

The effective pion decay constant /* is obtained by using the GT relation (5.17). Con- 
sidering the ratios f*/ fn (Table 5.5) and M^/Mn (Table 5.4) it may be concluded that they 
both decrease in the nuclear medium. It would be interesting to compare these ratios with 
well known Brown- Rho scaling law [101] f^/ fn = M^/Mi\f predicted by a simple Skyrme 
model. As it is seen from Tables 5.4 and 5.5 in our model f*//^ 7^ M^/Mn. Clearly, this 
deviation from Brown-Rho scaling law is caused by the the additional terms including the 
dilaton field. 

The masses of the omega and sigma mesons are supposed to decrease by the scaling 
law.There are experimental indications that this is true.In the same way the a — N coupling 
constant is expected to decrease. The ratio of g'*jvAr to its free value is presented in Table 
5.5. The changes in g^^NN and gaNN are nearly the same. Both are reduced in the medium 

by ~ 25% at p = po- 

In Figs. 5.11 and 5.12, the renormalized nNN and aNN vertex form factors respec- 
tively at p = po (dashed lines) in comparison with these in the free space (solid lines) are 
displayed. Appreciable quenching of both form factors is observed. At small momentum 
transfer these can be parametrized by a monopole form i.e. G-„NN{(i'^) = 9-kNN /{^ + 9*^/^^) 
and GaNNiq'^) — gaNN/{^ + 9 ^/A^). Table 5.5 shows that the cut off parameter A^ is 
decreased significantly at p = po- Relatively small changes in Go-ArAr(? ^) seem to be caused 
by a stiffness of the cr-field or equivalently C* [94] . 

The nucleon axial form factor (^^(^ ^)/Gyi(0)calculated for p = and p = po is presented 
in Fig. 5.13 with solid and dashed curves respectively. It is seen that, the modification of 
Ga{(i'^) is not as simple as that of G'7rArw(^^)- The medium leads to a quenching of the 
meson nucleon form factors over a range oi q^, while the quenching of Gyi(g ^) takes place 
at = and > 10/m~^. 

A further interesting quantity is the pion-nucleon sigma term E^riv- It is both the chiral 
symmetry breaking piece of the nucleon mass and a measure of the scalar density of quarks 
inside the nucleon. Due to the Hellman- Feynman theorem, it may be easily calculated [102] 
in the Skyrme model that 

For free space (p = 0) the Lagrangian including scalar mesons (5.8) gives S^at ^ 20.1MeV 

tViat is miipVi smallpr than that oKtainprl in thp nricrinal Skvrmp mnHpl fR.'^ R4l Tt mav Kp 
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Figure 5.11: The ttA^'A^ form factor. The sohd and dashed curves give results for the 
free space (p = 0) and the nuclear matter (p = po) respectively. 



argued that E^jv may also undergo changes in the nuclear medium i.e. E*^^ ^ E^jv- Actually, 
PCAC allows us to relate E^riv to the soft-pion limit of nN scattering [103] whose parameters 
may not be the same in free space and the medium. So, defining an effective value of the 
in-medium pion-nucleon sigma commutator as 



m *{N\ drt^\Ny 



*2 



(5.25) 



where in the used notations |A^)* - is the state of the nucleon bound in nuclear matter, 
we obtain E*^^ ^ AO.lMeV at normal nuclear density. This means a large increase of the 
nucleon sigma term: YT^^/Ti.^i^ 2. Using appropriate solutions of Eqs. (Al.13), it is 
estimated that 



m 



*2 



■'■kN , 



-"kN 



J e-^^'[l-cos{e*)ydx 



= 1.^ 
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ml J e-^^ll - cos{e)]x^dx 



mi 



(5.26) 
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Figure 5.12: The crA^'A'" form factor. Solid and dashed 
curves are for p = and p — Po respectively. 

the medium renormalization of J2t:n is caused mainly due to a large modification of the 
profile functions (see Figs. 5.14 and 5.15). On the other hand Ettat gives a good estimate for 
the quark condensate at finite density: 

vac T^-K fn 

Hence, assuming the last equation holds it is concluded that the in medium enhancement of 
EjrAT leads to further quenching of the scalar density of quarks {qq)p in nuclear matter. 

5.7 Discussions and summary 

We have proposed a medium modified Skyrme like Lagrangian which takes into account 
the distortion of the basic nonlinear meson fields by the nuclear medium. It has been 
extended by including the scalar - isoscalar sigma meson. The influence of the medium on 
pion fields is introduced by the self energy operators flp, (= Ha) and fig while the effect of 



(5.27) 
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Figure 5.13: The axial form factor. Solid and dashed curves are for p = and p = po 
respectively. 

The Lagrangian is applied to study changes in the hadron masses and meson - nucleon vertex 
form - factors. 

In particular, the mass of the A - resonance decreases more than that of the nucleon in 
the nuclear medium. Consequently the pion requires lesser energy to excite the nucleon to 
the A in the nuclear medium than it does for a free nucleon. The mass difference between 
A and N decreases to 42% of that for free particles at the nuclear matter density. This is 
quite consistent with earlier estimates [93, 104] but contradicts the recent theoretical results 
of Mukhopadhyay and Vento [105], who found SM^^/SM^^ pa 1.25 for p — O.Spo- So, it 
would be quite interesting to study NA mass splitting experimentally by an analyses of the 
NA transitions in heavy nuclei. 

We have investigated the characteristic changes of the decay width Fcr^TTTr at zero tem- 
perature. One may expect that the temperature-dependence of the physical quantities is 
qualitatively similar to the p dependence. In this sense, our results are in good agreement 
with predictions of the in - medium NJL model [106], that at sufficiently high temperatures 
the a meson becomes a sharp resonance and its width may even vanish. Clearly, more precise 
predictions on F^^^^fT, o) in the framework of the present Laerangian should be made bv 
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Figure 5.14: The profile function 9{r) for the p — (soUd curve) and p — po (dashed 
curve). They are the solutions of equations of motion (Al.13) in the sector with 
S = 1. 

studying thermal Green's function of a - meson e.g. within Thermo Field Dynamics. 

Furthermore the in - medium version of GT relation (5.17) holds in the present approach. 
The renormalized pion decay constant /* and nucleon mass do not satisfy the BR scaling 
[101]: the change in the nucleon mass is larger than that in /j^. 

The medium effects lead to a quenching of the meson - nucleon form factors as well as 
the coupling constants g^NN and QaNN- The latter should be compared with the results of 
Banerjii and Tjon [107] obtained in the framework of CCM in the meanfield approximation. 
There Qcnn increased at low densities making a looping behavior of the binding energy at 
saturation as a function of the density. However we believe that the natural reduction of the 
meson masses and coupling strengths found in the present model are expected to give a good 
description of the saturation properties of nuclear matter. So, it would be quite interesting to 
study the saturation properties of nuclear matter in relativistic Brueckner approach by using 
density dependent meson masses and coupling constants obtained in the present model. 

By introducing a formal definition of the in - medium pion - nucleon sigma term S*^ ~ 
m*^\dMN* /dmt^], it is found that in contrast to the meson - nucleon couphng strengths the 
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Figure 5.15: The profile function a{r) for the p = (sohd curve) and p — po (dashed 
curve). They are the solutions of equations of motion (Al.13) in the sector with 

Ettat increased in the medium: E*^/E^Ar ~ 2. This enhancement could lead to a decrease 
of the quark condensate {qq)p in nuclear matter. However as it was recently pointed out 
by Birse [108] this change should not lead to a drastic and rapid restoration of the chiral 
symmetry in nuclear matter. 

It is anticipated that ultrarelativistic heavy - ion collision experiments (e.g. at RHIC) 
will provide significant new information on the strong interactions through the detection of 
changes in hadronic properties [17, 109]. This would provide an impetus to consider refined 
models for strong interactions in nuclear matter. 

In conclusion, the modified in - medium Skyrme Lagrangian with scale invariance provides 
a useful insight into the role of medium in changing various properties of the mesons and 
nucleons. Even though scale invariance is badly broken in strong interactions its inclusion 
gives important information on the role of this symmetry property in a many-particle system. 



Chapter 6 



Nucleon electromagnetic form - factors in a nuclear 
medium at zero temperature. ^ 

6.1 Introduction 

The number of studies devoted to in - medium renormalization of the properties of nu- 
cleons has ever been increasing in recent years. Despite the absence of precise data from 
experiments from heavy - ion coUisions and deep inelastic nucleon scattering, there are rea- 
sons to believe that nucleons immersed in a nuclear medium undergo swelling. Theoretically 
it is clear that this internuclear effect must be due to the effect of strong meson fields on the 
quarks constituting a given nucleon and on its mesonic cloud. 

It is natural to assume that the change in the nucleon dimensions due to the intra 
nuclear meson field is of a more general character. Specifically nucleon elastic and inelastic 
form factors in a nuclear medium are expected to be different from those in a vacuum. 

An effective Skyrme Lagrangian describing the properties of a nucleon immersed in a 
nuclear medium has been given in the previous Chapter. Here we analyze the nucleon 
electromagnetic form factors in a nuclear medium. 

6.2 Electromagnetic properties of nucleons 

The problem of calculating of electromagnetic form factors for the nucleons reduces to 
computing isoscalar and isovector currents in a model. In topological models, in particular, 
in the Skyrme model, the isoscalar current is proportional to the baryonic current 

B^=^e^''''^TrL,L^L0, (6.1) 

which is seen to be model independent. This expression remains in force in the theory based 
on the Lagrangian in Eq. (5.8). Since the isoscalar mean square radius {r^)j^Q is expressed 
in terms of the zero component of the baryon current as 

oo 

{r')j=o = J B'd'r , (6.2) 



variations in it are not very important. Nontrivial modifications are due exclusively to the 
distortion of the skyrmion profile function in the nuclear matter. Similarly, the isovector 

^The Dresent chanter is based on following articles bv the author and his collaborators: [78, 79, 80, 811 



mean - square radius {r'^)j^i which is determined by the zero component of the vector 
current 



= -i^C^T^ AW + ^m) + r{[U[L^, L,]] + [R^R^, R,]]} , (6.3) 



where 



La = U+dM; Ra = UdM+; C.^l , ^ ' ^ ° 

^ ^ ^ I (1-Xa) , = 1,2,3 



[the vector current is treated as a Noether current that corresponds to the transformation 
U{x) exp{iQi^)U{x) cxp(— jQi?) {Ql,r -2x2 matrices], does not depend exphcitly on the 
medium parameters either. Relevant expressions for (r^) are presented in the Appendix B. 

Proceeding to analyze the isoscalar and isovector magnetic moments, we note that the 
former. 



does not depend on the medium parameters either. In contrast, the latter, /X7=i, which is 
expressed in terms of the spatial component of the vector current as 



involves the features of the nuclear medium explicitly. In our approach these stem from the 
contribution of the kinetic term in equation (5.8). 

The expressions for the electric and magnetic form factors, Ge and Gm respectively, 
can easily be obtained from equations (5.8) and (6.3) by using the technique developed in 
[18, 84]. The results are 




(6.4) 
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(6.7) 
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where 




(6.8) 



is the skyrmion moment of inertia, x = eF^r and jo) Ji are spherical Bessel functions. Note 
that, in nuclear matter, the form factors do not obey the well known relations [18, 84] 

Am oq 

(6.9) 

because the Lorentz invariance does not hold here. 



6.3 Input Parameters 

In our calculations we set the input parameters as = l%QMeV and Cg = (260MeV^)^, 
whereby the correct value is reproduced for the sigma meson mass ma — hhQMeV . The 
Skyrme parameter e was chosen in such a way as to fit the nucleon and delta - isobar masses 
in vacuum ( = 938MeK and Ma = 12?,2MeV respectively). The effect of nuclear 
medium was specified by fixing the parameters of the polarization operator Ha at g'^ — 0.6 
and Co = 0.13m~^ [85] as it was done in Chapter 5. 

Let us discuss the density dependence of the input parameters. As it was indicated in the 
previous chapter the renormalization of the Skyrme parameter e cannot be considered here 
because the p mesons have not been included on the effective Lagrangian (5.8) explicitly. 
For this reason we set e* = e. The renormalized pion mass m* has the form [86] 



. ;i - 0-22p/po 
yjl - 0.29p/po 

where po = 0.5m^ - is the normal density of nuclear matter. This parameterization also 
takes into account the correction to the contribution of the S - wave component of the pion 
self energy lis- 

The only input parameter in the dilaton sector is the gluon condensate C*, whose in - 
medium renormalization has to be clarified conclusively. In the proposed model, this quantity 
is given by Eq. (5.21) and can be expressed in terms of sigma meson mass as 

^9- (6.11) 

Unfortunately, no information about the renormalization of 171^ has been deduced so far 
within the present approach. In view of this, the parametrization 

m; = [1 - 0.12— (6.12) 
Po 

which was obtained in [89] within the quark - meson coupling model, is taken for the sigma 
meson mass. 



6.4 Results and Discussions 



Table 6.6: Ratio of the nucleon mass, nucleon root -mean - square radii, and nucleon mag- 
netic moments in a nuclear medium (labelled with an asterisk) to the corresponding vacuum 
values. 



p/po 











{rXn f^*I=0 f^*I=l P'l f^*n 


Mm \ 








(r2)^ lU^o 



0.0 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 

0.5 0.900 1.101 1.069 1.079 1.044 0.856 1.081 1.052 1.119 

1.0 0.819 1.196 1.122 1.145 1.061 0.817 1.145 1.103 1.200 



Table 6.6 presents the calculated root -mean - square radii and magnetic moments of the 
nucleons in the nuclear medium. It can be seen that, in - medium modifications to the root 
-mean - square radii of the proton are greater than those to the root -mean - square radii of 
the neutron. The effect is inverse for the magnetic moments - that is , the neutron magnetic 
moment is more sensitive to the impact of the medium than the proton magnetic moment. 

In nuclear matter of normal density, the proton charge radius (r^)^/^ increases by some 
14%, which is in agreement with the result presented in [111, 113]. We note that, only 5% 
increase in this quantity was found by Lu et al. [114], who took into account solely the 
renormalization of the nucleon mass, disregarding the renormalization of the pion - nucleon 
coupling constant g-n^NN- However, in previous Chapter it has been shown that, the in 
- medium modifications of meson - nucleon coupling constants are less sizeable than the 
corresponding modifications of the nucleon mass. 

Prom Table 6.6, it can be seen that the effect of a nuclear medium leads to an increase in 
the nucleon magnetic moment. As it was first revealed in [111, 112], the isoscalar magnetic 
moment changes more pronouncedly than the isovector one. In our calculations, the former 
and the latter change by about 14% and 8% respectively, at p = 0.5po and by 18% and 15% 
respectively, at p = pQ. 

Figures 6.16 and 6.17 display the proton and neutron electric form factors respectively, 
at the nuclear density values ofp = 0, p=0.5po and p — po (the solid, dashed and dotted 
curves respectively). 

It can be seen that with increasing the form factors decrease much more slowly in a 
nuclear medium than in a vacuum. Figures 6.18 and 6.19 show the normalized magnetic form 
factors G*jy,^{q'^)/G\,j{0) for the proton and neutron respectively. At ~ 0.3GeV^ (7.5-Fm^^) 
the proton and neutron charge form factors decrease by about 9% and 2% respectively at 
p = 0.5po and by about 17% and 4% at normal nuclear density. Similarly, the proton and 
neutron magnetic form factors decrease by about 6% and 5%, respectively, at p = 0.5po and 
by about 11% and 9% at normal nuclear density. 

From a comparison of the data in Figs. 6.16-6.19 it can be seen that the charge form 
factors are more sensitive to medium effects than the magnetic form factors. Qualitatively 
this conclusion complies with the results obtained in [114] on the basis of the quark - meson 
coupling model. 

We note , however, that at the above momentum transfer value that is at ~ 0.3GeV^ 
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Figure 6.16: Proton electric form factor The solid, dashed and dotted 

curves correspond to the nuclear - density value of p = 0, p = 0.5po and p = po 
respectively. 

Table 6.7: Form factors versus the nuclear density (p = Xpo) at some values of the momentum 
transfer Here the magnetic form factors are normalized to their values at q — 0. 
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0.087 


0.655 
0.618 
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0.650 
0.620 

0.590 
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Figure 6.17: The same as in Fig. 6.16 but for neutron electric form factor. 

p = 0.5po and by about 8% at normal nuclear density. A similar situation is observed for the 
magnetic form factors as well. That Lui et al.[114] found so weak suppression of the form 
factors in nuclear matter stems from their disregarding of the renormalization of meson - 
nucleon coupling constants and masses in their study. 

In ref. [115] the form factors for bound nucleons were calculated on the basis of the bag 
model. By preliminary fitting the input parameters of the model, the experimental values of 
the root - mean - square radii and magnetic moments were reproduced there rather well for 
the free nucleons. Further, medium effects were taken into account in that study through the 
renormalization of meson and nucleon masses. The analysis performed in [115] revealed that 
the nucleon electromagnetic form factors are suppressed in the ^^C, ^°Ca and ^^Fe nuclei. 
The nucleon magnetic moments in medium were enhanced by 2% to 20% for various nuclear 
species. These results agree with ours and the ones from refs. [Ill, 112] However, the proton 
charge radius in the ^"Ca nucleus is enhanced in ref. [115] by nearly 40%, which seems too 
large in comparison with our results (about ~ 15% at the normal nuclear density). It should 
also be noted that the results reported in [115] do not show the suppression of the axial 
coupling constant in nuclei. The failure to reproduce this well known effect casts serious 
doubts on the nonnent behind the calculations from f1 15l. 
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Figure 6.18: Normalized proton magnetic form factor G^((7^)/G'^(0). The solid, 
dashed and dotted curves correspond to the nuclear - density value of p = 0, p = 
0.5po and p — Po respectively. 



It is well known that, in a vacuum, the neutron and proton form factors calculated in 
topological models at low momentum transfers obey following relations [12, 18] 

Our calculations have demonstrated that, because of violation of Lorentz invariance, this 
relation is fulfilled only partly: 

We have also established that, at low momentum transfers {q^ < l.bFmr'^), the form factors 
show a nearly linear decrease with increasing a nuclear density. This results are in line with 
the results obtained in refs. [Ill, 112]. The nuclear density dependence of the form factors 
is illustrated in Tabl. 6.7 for various values of the momentum transfer. 
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Figure 6.19: Normalized neutron magnetic form factor G^(g^)/G^(0). The solid, 
dashed and dotted curves correspond to the nuclear - density value of p = 0, p = 
0.5po and p — Po respectively. 

6.5 Conclusions and Summary 

On the basis of chiral and scale invariant Lagrangian, we have investigated the possible 
modifications to the nucleon electromagnetic form factors in a nuclear medium. Our results 
demonstrate that this model faithfully reproduces the features of hadrons as functions of the 
medium density. The calculations have revealed that the isoscalar nucleon form factors are 
more sensitive to the medium effects than the isovector ones. On the whole, our predictions 
are consistent with the results of other authors. This proves that, qualitatively, the resulting 
behavior of the features of nucleons in a nuclear medium is model independent. Specifically 
the coupling constants, masses and form factors are suppressed by the effective fields of the 
nuclear medium, wchcrcas the dimensions and magnetic moments arc enhanced. By way 
of example we indicate that, at the normal nuclear density the root mean square radii for 
protons and neutrons increase by some ~ 15% ~ 6% respectively. With increasing density 
the nucleon form factors exhibit a nearly linear decrease in the region < IFm"^. 

Here we have taken no account of recoil effects. Both prior to and after coUision with 



a photon, the nucleoli involved is assumed here to be frozen at the center of nucleus where 
the density is constant. As it was shown by Bunatian [116] the contribution to {r'^)N from 
recoil effects is small, but it can be expected that modifications to the nucleon structure will 
greatly depend on the nucleon momentum within Fermi sphere. The inclusion of recoil effects 
in the calculation of the nucleon electromagnetic form factors for finite nuclei would be of 
great interest in the problems of electron nucleon scattering. Calculations of this type my 
also be helpful in processing data from future MIT Bates and TJNAF experiments [117] that 
will measure the neutron electromagnetic form factor. Since there is no neutron free target 
relevant information will be extracted from data on electron nucleus (say e ^He) scattering 
where it is necessary to take onto account nuclear medium effects. 

Note that, the modification of nucleon electromagnetic form- factors in finite nuclei has 
been considered recently in ref. [118]. 



Chapter 7 



Meson - Nucleon Vertex Form Factors at Finite 

Temperature. ^ 

7.1 Introduction 

In-medium properties of hadrons and their interactions is a field of high current inter- 
est. For studying phase transitions and thermodynamical properties of the nuclear system 
under extreme conditions, it is essential to determine the temperature (T) and density (p) 
dependence of the nuclear force. Below a critical temperature Tc, meson exchange picture of 
nuclear physics provides the natural description for the nucleon - nucleon (NN) interaction. 
Quarks and gluons are certainly present in all A^A^ interactions, but it is not always neces- 
sary to take them explicitly into account especially below a critical temperature T^, when the 
transition into a quark gluon phase takes place. At distances smaller than a fermi, the inner 
structure of hadrons is probed since it involves the short distance or the high momentum 
components of the wave function. This structure is typically taken into account by including 
vertex form factors. In hot and dense matter, the structure of hadrons undergoes changes 
which should lead to a modification of the meson - nucleon coupling constants or, in general, 
the form factors. The temperature [121] and density dependence of some coupling constants 
have been investigated recently (see Chapters 5, 6), while that of form factors is still fairly 
unknown and will be studied in the present Chapter for the first time. For this purpose we 
shall use the formahsm of Thermo Field Dynamics (TFD) in order to obtain the required 
temperature and density dependence. 

The Thermo Field Dynamics (TFD) which was first suggested by Takahashi and Umezawa 
[122] - [126] is a real time operator formalism of finite temperature quantum field theory. In 
this framework, all the operator formalism of quantum field theory at zero temperature can 
be extended directly to finite temperature and density. Therefore in TFD it is possible to 
continue using Wick's theorem and the Feynman diagrammatic approach as in the case of 
zero temperature theory. An auxiliary field is introduced which leads to a field doublet; the 
propagators become 2x2 matrices, the Feynman rules are now algebraic operational rules 
in the space of 2 x 2 matrices. Recently, some attempts to study in medium properties of 
hadrons at finite temperature using TFD [127, 128] have been made. 

A study of relativistic heavy-ion collisions [129] is expected to lead us to a proper equation 
of state of hot and dense nuclear matter. The behaviour of nuclear matter as a function 
of temperature and density is relevant for investigation in nuclear physics, astrophysics, 

^The present chapter is based on following articles by the author and his collaborators: [119, 120] 



cosmology and particle physics. The phase transition from hadronic system to a quark- 
gluon plasma phase and the subsequent hadronisation of the quark-gluon plasma phase can 
provide information on the asymptotic freedom as well as confinement behaviour in Quantum 
Chromodynamics (QCD). Therefore it is important to study theoretically the behaviour of 
the hadronic system at high temperatures. In this Chapter we shall make an attempt to 
study such a behaviour. In Sec.2, a brief review of the TFD is given. In Sec. 3 we calculate 
the renormalization of meson nucleon vertices using the Feynman propagators from TFD. 
The results of calculations and their discussions will be presented in Sec. 4. 

7.2 Finite temperature formalism for meson and nu- 
cleon propagators. 

The Thermo Field Dynamics is a real time operator formalism of quantum field theory 
at finite temperature. The main feature of the TFD is that thermal average of operator A is 
defined as the expectation value with respect to a temperature dependent vacuum, \0{f3)), 
which is obtained from the regular vacuum by a Bogoliubov transformation. Therefore, we 
have 

< A >= Tr(ie-^(^-'^))/Tr(e-^(^-'^)) = (0(/3)|i|0(/3)) , (7.1) 

where P = l/ksT with ks being the Boltzmann constant, H is the total Hamiltonian of the 
system, and fj, is the chemical potential. Due to the doubling of all degrees of freedom, every 
field is represented in the thermal doublet form in TFD. For example, the thermal doublet 
for the nucleon field is 

^'»'(-) - { ] . (7-2) 

where iIj{x) is the ordinary nucleon field and ip is the doublet partner of iIj{x). Here the 
superscript represents the transpose operation on the vector index. As a consequence, the 
thermal propagator of a fermion field is a 2 x 2 matrix defined as 

zSp'\x„x,) =< 0(/5)|T[^(»)(xi)V'('')C^2)]|0(/?) > (7.3) 

The free thermal propagator Sp to be used in perturbation theory can be separated 
into two parts, i.e. the usual Feynman part 5*0, propagator at zero temperature, and the 
temperature dependent part St such that Sp''^ = Sq"'''^ + S^'^K For a fermion with mass M, 
the propagators have the form 

.r = (MM),.,(«f' (7.4) 
4"" - 2.iS(,^ - M^)i^+MU ( iX?^; ) (7.5) 



with 



G'o(p^) = ^ I • ' 

— + 16 



(7.6) 



where x = P{po — lJ'), a and b are Dirac indices, and 9{x) is the step function. The propagator 
of the scalar particle with mass is given by 

A^(p)=Ao^(p) + Af(p), (7.7) 

(7.9) 

Dr = -2ni5{p^ - m|) ( f "^Jf,- ^^^"f>> ] (7.10) 

y isinh20po sinhVpo / 



with 



i^o(p^) 



1 



— + ie 



= (i_e-M)i/2 ' (7-11) 



sinh0, 



,-|2/|/2 



(l_e-M)V2' 

where y = Ppq. Note that the propagator of vector mesons would be similar to that of the 
scalar mesons except for an additional factor of {—g^iv + k^ky/rn?). 

7.3 Vertex form factors 

The OBE model of the N-N interaction [62] includes exchange of several mesons. Here 
we consider the meson-nucleon interactions of tt , o" , a;- and p-mesons: 

^pNN = [-gpNNi^lu.T'il^P^ - ^^^iJTaf,^^{d^p^ - d^Pf,)] . 

Here gsNN are the meson-nucleon coupling constants and Jojnn and fpNN are the tensor 
coupling constants for cu- and p- mesons respectively. 

The corresponding meson nucleon form factors are usually defined as: 

{N{p')\T^\N{p)) = ig^NN{t)u{p')^5rau{p) 
{Nip')\T^\N{p)) = -g.NN{t)u{p')u{p) 

{N{p')\V^jN{p))=u{p') 
(iV(p')|r"''^|iV(p)) = «(p')r„ 



u{p) (7-13) 
u{p) , 



where q-^nni QuNNi G^nn, (F'^nn) , Gp^N and {Fpj^N) are form factors of meson - nucleon 
vertices; t = = — [qf = (p — p')^ is the 4-momentum transfer and M is the nucleon 
mass. 



p' p' p' 




Figure 7.20: Fcynman diagrams for pion - nucleon vertex. The solid line is for 
nucleon. Dashed, dotted, dot - dashed and wavy lines are for tt, ct, uj and p mesons 
respectively. 



In order to investigate the T dependence we calculate three line vertex correction, as is 
illustrated in Fig. 7.20 e.g. for gT^i^i^{t,T). In accordance with Feynman rules this may be 
written as: 

TA{t,T)=TA(t)+lY,^AB{t,T), 

B 

(7.14) 

where A,B - denote n.a.uj.p mesons, e.g. F^ = gnNNit), a = p' — k.b = p — k. Here 
5'i?(a) and A^(A;^) arc the in - medium nucleon and meson propagators respectively. In 
the framework of TFD the thermal propagator has a 2 x 2 matrix structure, but only the 
11-component refers to the physical field. So, we may use the following representation: 

S}r\a) = Spia) = (a + M)[Go{a) + Gria)] , 

^o{k')^-r-2 DT{k') = -2mnB{k)6{e-m'), ^^'^^^ 

= r^^— - , Gria) = 27ri5{a^ - M^)NF{a) , 

— -\-i£ 

where a = 0^7'^, Np[a) = 9{ao)nF{a) + 9{—ao)nF{a), m is the mass of corresponding meson 



and ^ ^ 

^^(«) = e/3(|ao|-M) + 1 ' = e/3(|ao|+/.) + 1 

(7.16) 

are the fermion, antifermion and boson distribution functions respectively. It is clear that 
when one substitutes Eq. (7.15) into Eq. (7.14), AAsityT) will be separated into two parts: 
one refers to the naive zero temperature contribution, while the other one depends on the 
density and temperature [127]. We shall concentrate on the latter part, which may be 
rewritten as follows: 

AAB{t,T)= f ^TB{k'){d + M)TAmb + M)TB{k')x 

^ ^2^) (7.17) 

x[Go{a)Go{b)D^{k^) + 2Go{a)GT{h)D^{k^)] 

Note that wc have neglected terms that are quadratic in the temperature dependent distri- 
bution function. Now, using (7.17) and (7.15) in (7.14) we get the following expressions for 
the form factors: 

GANN{t,T)/GANN{t,T = Q) = l + lY, [ ^^WAB{t,k\T)x 

° ^ ^'^^ (7.18) 



B 

B/u2 



x[Go{a)Go{b)D^{k^) + 2Go{a)GTib)D^{k 

where the explicit formulas for WAB{t,k'^,T) may be found in the Appendix C. Here we 
present the expressions for Wui-n-it, k^, T) and WpT^{t, k^, T) for illustration: 

WUt,k\T) = |2[2M^ - {ab)] - l^^} , 

(7.19) 

W^it^^T) = Z^k^ j2[2M^ _ (a6)] - j , 

where W^T^{t, k"^, T) and WpT^{t, k^, T) denote the contribution from the tt- exchange diagram 
to the Gu,NN{t,T) and GpNN{t,T) respectively. 



7.4 Results and discussions 

In order to start the calculations, a set of free space meson - nucleon form factors are 

chosen. We choose the OBE monopole form factors [62] (Bonn A): GsNNit) — gBNN{^\ — 
m|)/(A| - t) ,where, for example, gl^j^/An = 14.09, A^nn = 1005Mey, m„ = 550MeV. 
Let us first discuss the T dependence of meson - nucleon coupling constants - gBNN{T) = 
GsNNit = nT-^jT). The variation of the coupling constants g-KNN{T) / g-KNN{T = 0), 
gaNN{T)/gaNN{T = Q)-, gu:NN{T) / g^NN{T = 0) and gpNN{T) / gpNN{T = 0) with temperature 
are displayed in Fig. 7.21. It is clear that they are nearly independent of the tempera- 
ture below T^, then change rapidly for T > . A similar behavior was also predicted 
by Zhang et al. [127] and by Dominguez et al. [121] for pion - nucleon coupling constant 




Figure 7.21: The ratio of meson-nucleon coupling constants at finite temperature T 
at T=0 as a function of temperature; p — 0. 



g-KNN{T)/ g-^NN{T — 0). But here in Fig. 7.21, an unexpected result is that, the uj — N cou- 
phng constant, gojNN{T) increases while the coupling constants of all other mesons decrease! 
Let's consider a possible origin of this controversy, by comparing the medium modification 
of uNN and pNN coupling constants. In the present model the modifications arise from the 
triangle diagrams (Fig. 7.20). Actual calculations show that, for any GBNN{t,T) the trian- 
gle diagram with pion exchange gives a dominant contribution especially at small density. 
It is clear from Eq. (7.19) that at small t the explicit expressions for gujNN{T) and Qp^NiT) 
formally coincide. However, Wp^, = and hence the two have an opposite sign. The 

factor (—1/3) arises from the isotopic spin structure. In fact, using Eq. (7.14) and Fig. 7. 20, 
gu!NN{T) and gpNN{T) may be written in a schematic way: 

gujNN{T) g^NN{T = 0) + fif^jvAT XI Tpg^jMNTp + . . . 

P (7.20) 

gpNN{T)Ta ^ gpNN{T = 0)t„ + g'^jviV X] TpTagpNNTp + ... 



where a and (3 are isospin indices. 




T (MeV) 



Figure 7.22: The ratio of mean square radius at finite temperature T and at T = 
as a function of temperature; p — 0. 

We omit the spin variables, since uj and p have the same spin structure. Now, using y^r^ = 3 

and ^^TpTaTp = —Toil it is clear that the contribution of the leading triangle diagram with 

P 

pion exchange to guiNN{T) and gpNN{T) has different signs. In other words, the term A^^r in 
Eq. (7.14) is negative, while K^-,, is positive. 

In general, we conclude from the results (Fig. 7.21) that at a certain critical temper- 
ature , the couplings Qbnn change drastically. The thermal behavior of gnNNiT) = 
gnNN{t,T)\t-fO has been investigated earlier [121] and we get here a similar behaviour of 
g-KNN{T) as a function of temperature. The temperature Tc, where g^NNiT) changes dra- 
matically, was interpreted in [121] as a signal for the quark - gluon deconfinement phase tran- 
sition. Indeed, near Tc the associated mean square radius < r^jy^y >= 6[^ In GsNNit, T)][t_>o 
is a monotonically increasing function of T and in fact diverges at the critical temperature. 
A similar behavior of < r^jvAr > i^ ^^^^ found in our calculations and is illustrated in Fig. 
7.22. As the critical temperature, T^, is approached the strength of coupling of vr, a and 
p mesons to nucleons is quenched, at the same time, the size of nucleons as probed by the 



Tabic 7.8: Parameters of vertex form factors in Eq.s (7.21)-(7.22) at half of normal nuclear 
matter density p — O.Spo- 
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Table 7.9: The 
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mesons: T^ ^ 360Mey, T^ ^ 95MeV, T^ 175Mey and TP ^ 200Mey. In OBE picture 
this means that , in the temperature region e.g. 200Mey < T < 300Mey, the p meson 
exchange is no longer important while the pion exchange is still important. On the other 
hand, the a and vr mesons are mainly responsible for the attraction between nucleons. So, 
the quenching of the aNN coupling constant at T > T^ leads to the vanishing of the bound 
state, as it was predicted earlier [127]. 

The density dependence of vertex form factors has been studied in greater detail [81, 93]. 
It was found that most of the vertex form factors are quenched at high densities and it was 
anticipated that the temperature dependence is likely to yield results that are qualitatively 
similar to those of density dependence. 

Now, let's consider their temperature dependence in some detail. The form factors as a 
function of momentum transfer at several temperatures are displayed in Figs. 7.23 - 7.25. 
It is seen that, the in-medium effects lead to the suppression of gTTNN{t,T), QaNNitjT) and 
GpNN{t,T). The temperature dependence of Gi^i\fN{t,T), is opposite to that of Gp^^ due to 
it's isotopic structure as outlined above. The temperature dependence of tensor couplings 
of vector mesons FvNN{t,T) are quite similar to that of vector couplings GvNN{t,T). Par- 
ticularly, the ratio Ky — Fvnn/Gvnn, where Kp — 6.1 and = in free space, remains 
constant in a wide range of temperature. For practical calculations a parametrization of these 
GBNN{t, T, p) form factors is needed. At small momentum transfer we can parametrize them 
by a monopole form : 

GBNN{t, T, p) = gsiT, p){AUT, p) - ml)/{Al{T, p) + t) (7.21) 



Table 7.10: The same as in Table 7.8 , but for p — 3po 
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Figure 7.23: Meson nucleoli form factors at several temperatures as a function of (f 
using parametrisation of BONN group [62] for ttA^A^ vertex at p = 0. 



where in general the effective mass of a meson, m^, is also temperature dependent. Consider- 
ation of this dependence is beyond the scope of the present Chapter. But here, for simplicity, 
one may consider this as just a parameterization and choose the parameters gsiT, p) and 
A^(T, p). Their temperature dependence is still unknown. Here the T dependence may be 
represented, for T << T^, in a polynomial form as: 



9BiT,p) 



^b(T = 0,p = 0) 

Ab{T,p) 
As(T = 0,p = 0) 



(7.22) 



where $(p) = 1/(1 + Cop/po)- The calculated form factors arc fitted to this form and the 
paramters Co, a and /? arc determined. The results are presented in Tables 7.8, 7.9 and 7.10 
for densities p = 0.5pQ, p = po and p = 3po, respectively (with po = 0.17/m^'^ - the density 
of normal nuclear matter). The results show that the density dependence of parameters, 
and hence, the form factors is not so drastic fsharnl as their temoerature deoendence in 
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Figure 7.24: The same as in Fig 7.23 but for aNN vertex. 

the range of moderate densities. So, at T the Eq. (7.22) with Cq = 0.26 is in good 
agreement with Brown Rho scahng law [131] (where Co = 0.28). 

Summarizing this Chapter, we have considered the temperature dependence of meson 
- nucleon form factors and couphng constants in TFD formahsm. It is shown that at a 
critical temperature, where the coupling constant changes drastically and the associated 
mean square radius diverges is different for different mesons. The temperature dependent 
vertex form factors are parametrized in a simple monopole form and the T-dependence of 
these parameters is clarified. These form factors may be used in the calculation of the in - 
medium NN cross sections [130] and in investigations of the properties of hot dense matter. 

In the present calculations, the monopole form factors with the parameters (at T = 0) 
given by the Bonn group [62] are used. But the next question which may arise is: are the 
results sensitive to the shape of the input form factors? To answer this question a different 
set of form factors obtained in Chapter IV in the framework of a topological soliton model 
have been used. The thermal behaviour of such form factors is shown in Figs. 7.27 - 7.30. 
Now, by comparing Figs. 7.23 - 7.26 and Figs. 7. 27 - 7.30 it is clear that only the aNN form 
factor is affected the most. The reason is that even in free space and at T = the aNN 
form factor in the tonolomnal soliton model f54l is munh harrier fA_ Ri (^OOMp.V^ than those 
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Figure 7.25: The same as in Fig 7.23 but for uNN vertex. 



of the Bonn potential (A^ ~ 2000Mey) . Nevertheless, our main conclusions about thermal 
behavior of different meson nucleon vertices remain valid. 

The variation of the meson-nucleon coupling constant with temperature (Fig. 7.21) 
clearly indicates that the attractive NN force due to cr-exchange decreases quite rapidly 
while the repulsive N-N force due to cu-cxchanges increases. This would lead to the fact that 
the nuclear matter is quite likely un-bound at high temperature. In fact it will look more 
like a hard sphere gas since the repulsive interaction at short distances will dominate. The 
attractive interactions in such a case play only a small perturbative role. The hot and dense 
nuclear matter could be approximated as a free gas with an excluded volume around each 
nucleon. Perhaps the extensive studies [132]-[135] along this line can be justified on the basis 
of the results obtained in this Chapter. The nuclear matter would be hard to compress. 

It would be quite interesting to know, if the critical temperature, T^, maybe uniquely 
considered as a point of transition from the hadron state to a quark gluon state, as it was 
suggested [121]. The usual way of estimating Tc is based on QCD lattice calculations [136] 
or temperature dependent quark - gluon potentials [137]. In such calculations it is usual 
to start from high temperatures and then decreases the temperature, looking for a phase 
transition from the nuark Fluon state to a, hadron state i.e. hadronization of the nuarks. 
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Figure 7.26: The same as in Fig 7.23 but for pNN vertex. 

Clearly this method gives a unique Tc, since at any stage there is either a hadron or a 
quark gluon plasma phase. Contrary to such studies the present calculations explore the 
transition from a hadron to a quark-gluon phase and find that the critical temperature for 
hadron quark-gluon transition depends on the kind of hadrons under consideration. A 
similar result, where hadrons and quark - gluon plasma coexist have been obtained in lattice 
calculations [136]. In conclusion it is important to emphasize that the dynamics of hadrons 
and of quark-gluon plama at finite temperature is poorly understood. The question of critical 
phenomenon, in particular critical temperature, in the hadronic systems is not well-known. 
Studies in nuclear matter are needed urgently to clarify the critical transition from a hadronic 
system to the quark gluon phase which will depend on both density and temperature. The 
experiments with relativistic heavy ion colliders would need such an understanding to clarify 
the production of quark-gluon plasma from hadrons and the eventual hadronisation of the 
quark-gluon system. It is possible that there will be clear and unambiguous signals for the 
formation and subsequent hadronisation of the quark-gluon plasma [138, 139] . 
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Figure 7.27: Meson- nucleoli form factors at several temperatures as a function of 
using parametrisation of Meissnar et al. [54] for t^NN vertex at p = 0. 

G*t^nn{q^) G^jvAr(9^) used in Chapter 5. The small fluctuations around the vacuum 
value are related to the pion field by 

U = exp(2if7f/F^) fti 1 + 2iTT^jF^ - 2TejFl + . . . (Al.l) 

which gives the following approximation for the Lagrangian in Eq. (5.8): 

L^-\{SJi:fa,-\wr:^^ (Al.2) 

and for the equation of motion: 

-(VV)Q;p + m;2f = (A1.3) 

The factor ap = 1 — xa in the last equations is obtained by renormalization of the pion 
propagator in the medium i.e. by the A - hole self energy term Ha- We define the in- 
medium ttNN coupling constant and the vertex form factor by introducing the source term 
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Figure 7.28: The same as in Fig. 7.27 but for aNN vertex. 



into the equation of motion: 

-V^Tfap + m*^7r = iG*jvjvV^75f'0, (Al.5) 

where 

piPx / TP* I Tur* \ 1/2 / 1 \ 

where E*'^ = + M*^. In the Brcit frame {P' = —P, q = P — P') the matrix element of 
the source evaluated between nucleon states is given by: 



Using the quantization rules 
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Figure 7.29: The same as in Fig. 7.27 but for uoNN vertex. 



we have 



{N'{P')\7r^{f,R)\N{P)) ^ J {N'\7r^{x)\N)e-'^^^dx 



(A1.9) 



tq r 



( Va (rx) I N)e-'^^ sm{9) dx 



6(27r)3 

where 9 is defined by the hedgehog ansatz C/o(^) — e*^^'^^^^). Now the matrix element of 
Eq. (Al.5) is evaluated between collective wave functions for spin-up proton \p f) with 
momentum P and P' in the Breit frame and using equations (Al.6) - (Al.9) to obtain 
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(ALIO) 



3 J qx 
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Figure 7.30: The same as in Fig. 7.27 but for pNN vertex. 

with c = cos(^), s — sin(^). 

Similarly, the couphng constant at the aNN vertex is defined by the equation: 
(— + m*^)a = G^.pfj^'ilj'ijj. Evaluating matrix elements of both sides of this equation it is 
easy to obtain the aNN form factor: 

oo 

G^nAq^) = '^^F^J 3o{qx)S,{x)dx , (Al.12) 

o 

with Sa{x) = -x^a" + 2xa' + x^mfa. 

Note that the profile functions 9{r) and a{r) in St^{x) and Scr{x) are the solutions of the 
equations of motion: 

e"x^X^Oip + 4s2^'^ + 8s^e" + 2xe'x^ap + 2x'^e'xx'oip- 

-X^apS2 - 4s2d - x'^p\h = f A 1 1 

x\" + 2xx'-2xapx'{e'y2 + d)-16x'V,ff{x'-X'-')- ^ ' ' 

-3x2/?2(i - c)x^ = 0. 

where x = eF^r, S2 = sm{29), d = s^/x^, (5 = m^/eF^, Pg// = C*j2Ae^F^, 
X = exp(— ^(a;)). The boundary conditions are: 



/ r\ C r\ 1/1 / 1 , n \ — Rnf I 9 



For completeness we write also the explicit expressions for I = Iq + I^t- 

27r 



Io = ^js'{e-'^ + A{e'^ + d)}x'dx, 

"° oc (A1.14) 
^ [{e'^c' + 2d}x^dx. 



^rt — IT- 





Appendix 2 

Here some auxiliary formulas for Chapter 6 are presented. The isoscalar and isovector 
mean square radii are given by 

oo 

{r^)i=o = -2 / r^O's'^dr , 



oo 

lr'sV + ^e" + d)]dr (A2.1) 



I rV[x' + 4(e'2 + d)]dr 



where the notations 9, s etc are given in Appendix A. The nucleon mean square radii {r^)pn 
are expressed in terms of (r^)j^Q and (r^)^^^ as: 

We note that the expression for mean square radius in a medium is similar to that in a 
vacuum, a solution to the equation (A1.13) with medium parameters must be taken in 
former case. 

The proton and neutron magnetic moments are given by 

f^P,n — 2 [A*p,n + A*p,n ] ) 

CO 



where /j^b is the Bohr magneton. For the nucleon form factors we have 



E — 2^^E + ^E) ' - 'jy^E ~ ^e) > 



(A2.3) 



Appendix 3 



Here the explicit expressions for WAsit, k'^, T) introduced in Chapter 7 are given. 

W,, = G,{k^)[M^ - {ah)] , W,, = W,,\g^^g. , 

W^^ = A^[M2 - (ab)] + 3F^{k^)F^{k^)z , (A3.1) 



mf. 



4M2 



(A3.2) 
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2Ff{t) 
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^[^[(a6) + SM^] - 2{ak){bq) - 2{aq){bk)]^ ; 



(A3.3) 



r-=[Fr(A;2)]2[2-— ] + 
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2M^ 



r- = F-{k^)F-{k^) , 



3W 



W, 



pTT 



Wp, = -w, 



pTV \ G-n — >Ga 



^pw ^ pit; ) ^^PP ^^u}uj\u)—>p J 



(A3.4) 



i i 

where = G^NN{k'^), G^ = G^NNik"^), F^ = G^nn, F^ = F^nn , z = (kq) and = t. 
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